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INTRODUCTION 


PROBLEM  IDENTIFICATION 


Fuels  management  is  becoming  an  ever 
increasingly  important  and  more  diffi- 
cult task  to  accomplish  each  year. 
Frequently,  the  treatment  of  slash!/ 
in  partial  cuts  is  more  difficult  than 
treatment  in  clearcuts,  primarily 
because  consideration  must  be  given  to 
the  residual  stand.  Partial  cutting  is 
an  ambiguous  term  which  includes  any 
type  of  commercial  harvesting  other  than 
clearcutting,  ranging  in  degree  of  cut- 
ting from  a seed  tree  cut  to  a very 
light  salvage,  sanitation  or  selection 
cut. 

With  recent  public  criticism  of  clear- 
cutting,  various  types  of  partial  cut- 
ting are  receiving  much  more  widespread 
use.  In  addition,  as  land  managers  move 
toward  more  intensive  management,  par- 
tial cutting  becomes  more  widely  prac- 
ticed in  the  form  of  intermediate  cuts. 

The  Northern  Region  (Region  1)  can  be  used 
to  illustrate  the  increased  use  of  par- 
tial cutting.  Three  years  ago  partial 
cuts  composed  approximately  45  percent 
of  the  Region's  total  planned  sell  as 
compared  to  approximately  80  percent  in 
fiscal  year  1974. 

Recently,  the  importance  of  adequate 
slash  treatment  has  been  reemphasized 
and  the  use  of  extra  fire  protection  as 
a slash  treatment  alternative  in  partial 
cuts  has  been  deemphasized.  It  is  esti- 
mated that  in  Region  1 there  are  over 
400,000  acres  of  untreated  slash  in  par- 
tial cuts.!/  Admittedly,  part  of  this 
acreage  may  not  need  any  special  treat- 
ment; however,  from  the  standpoint  of 
fire  hazard  alone,  treatment  is  probably 
needed  on  a significant  portion  of  this 
total  acreage. 

The  Missoula  Equipment  Development  Center 
(MEDC)  began  a project  in  FY  '74  investi- 
gating fuel  treatment  systems  for 
partially  cut  stands.  The  following  is 
an  interim  report  which  describes  the 
Center's  plans  for  the  project  and 
results  of  the  first  year's  efforts. 


Jy  Slash  or  fuel  in  this  report  refers  to 
residues  such  as  foliage,  twigs,  branch  wood, 
bark,  rotten  wood,  cull  logs,  and  trees  unmer- 
chantable because  of  defect  or  size.  This 
includes  not  only  material  resulting  from 
harvesting  but  also  material  existing  prior  to 
logging. 


Initially,  a considerable  amount  of  time 
and  effort  was  expended  by  a 5-man  team, 
comprised  of  engineers  and  foresters, 
identifying  current  problems  in  fuels 
management  in  partial  cuts  and  deter- 
mining what  contribution  MEDC  could  make 
in  finding  solutions  to  the  problems.  A 
literature  search  was  conducted.  Field 
personnel  were  contacted  and  queried  on 
the  subject.  The  project  was  discussed 
with  personnel  involved  in  slash  projects 
at  San  Dimas  Equipment  Development  Center, 
the  Northern  Forest  Fire  Laboratory 
(Intermountain  Forest  and  Range  Experi- 
ment Station),  the  Forest  Residues 
Reduction  Project  (Pacific  Northwest 
Forest  and  Range  Experiment  Station), 
and  the  Coram  Logging  Study  (Intermountain 
Forest  and  Range  Experiment  Station). 
Efforts  have  been  made  to  assure  coordi- 
nation with  these  research  projects  and 
to  avoid  duplication  of  work.  In  addi- 
tion, the  team  visited  several  partial 
cut  sales  with  local  field  personnel  and 
discussed  their  slash  problems  on  the 
Lolo,  Kootenai,  Flathead,  and  Clearwater 
National  Forests  in  Montana  and  Idaho. 

It  was  determined  that  few  land  managers 
are  able  to  visualize  or  predict  the 
resulting  slash  on  a partial  cut  sale 
prior  to  harvesting.  As  a consequence, 
in  many  instances,  the  resulting  slash 
is  accepted  as  a gigantic  "surprise 
package"  after  logging,  and  satisfactory 
treatment  may  be  impossible.  On  the 
other  hand,  treatment  may  be  possible 
but  funding  may  be  insufficient  to  accom- 
plish the  task  if  an  adequate  allowance 
was  not  made  during  presale  planning.  In 
general,  adequate  data  and  planning  tools 
relative  to  fuel  treatment  systems  are 
not  available  to  help  accomplish  good 
presale  planning.  Since  slash  treatment 
on  Forest  Service  sales  is  considered  a 
logging  cost,  an  estimate  of  the  total 
treatment  cost  is  necessary  when  the 
sale  appraisal  is  prepared.  If  the  esti- 
mate is  too  low,  additional  funding  may 
be  sought  from  appropriated  monies  which 
are  extremely  difficult  to  obtain,  if  not 


2/  Based  on  a conversation  with  Jack  Alley, 
Div.  of  Timber  Mgt.,  Region  1,  US  DA  Forest 
Service. 
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impossible.  Even  if  appropriated  monies 
are  obtained,  they  may  be  inadequate  to 
accomplish  the  necessary  work.  Therefore, 
adequate  treatment  may  not  be  achieved. 
Ideally,  a land  manager  should  be  able  to 
predict  his  fuel  treatment  needs  prior  to 
harvesting  and  be  aware  of  feasible  treat- 
ment alternatives  and  their  associated 
costs . 

In  summary,  the  problem  is  that  most  land 
managers  have  no  reliable  or  accurate 
method  of  predicting  fuel  residues  .resul t- 
ing  from  a particular  partial  cut  sale. 
Further,  it  is  extremely  difficult  to 
estimate  the  type  of  treatment  that  is 
necessary  or  suitable  to  reduce  fuels  to 
an  "acceptable"  level;  a level  that  is 
generally  not  known  in  quantitative  terms. 
Also,  it  is  difficult  to  estimate  the 
associated  costs  of  a particular  fuel 
treatment  for  a particular  set  of  sale 
conditions.  Experienced  or  historical 
costs,  even  for  the  most  commonly  used 
slash  treatments,  are  reported  on  a per 
acre  basis  and  provide  little  reference 
to  variables  such  as  slope,  fuel  loading, 
residual  spacing,  etc.  As  a result,  the 
range  of  experienced  costs  for  a partic- 
ular fuel  treatment  is  extremely  broad 
and  the  average  cost  is  of  questionable 
value  for  planning.  Finally,  land 
managers  may  not  be  familiar  with  all 
treatment  methods  that  are  feasible  for 
a particular  sale. 


GOALS  AND  OBJECTIVES 

Based  upon  our  problem  analysis,  the  goal 
of  the  project  is  to  develop  a decision- 
making model  or  technique,  which  is 
better  than  currently  used  methods,  to 
predict  fuel  loadings,  fuel  treatment 
needs,  treatment  alternatives  and  associ- 
ated costs  prior  to  harvesting  for  a 
specific  set  of  sale  conditions.  Plans 
are  to  develop  an  operational  decision- 
making model  within  another  year. 

Eventually,  increased  utilization  of 
logging  wastes  may  eliminate  the  problem 
of  slash  disposal.  Considerable  work  is 
underway  by  the  Forest  Products 
Laboratory  and  others  to  find  ways  of 
increasing  utilization  and  developing 


markets  for  logging  residues.  But 
complete  utilization  is  not  yet  a 
reality  in  many  areas,  and  slash  dis- 
posal is  still  a problem.  So  we  plan  to 
provide  aids  to  the  fieldman  that  can  be 
used  in  making  decisions  to  manage  slash 
in  partial  cuts  as  long  as  the  problem 
exists. 

Because  useful  fuel  information  from  past 
studies  is  available  for  tree  species  in 
Region  1,  and  it  is  easier  financially 
and  more  expedient  to  contact  field 
personnel  in  this  locale,  the  study  is 
initially  being  limited  to  the  Northern 
Region.  However,  we  expect  that  many 
techniques  developed  specifically  for 
this  area  may  be  applicable  in  other 
Regions  with  little  or  no  modification. 

Specific  objectives  which  will  lead 
toward  the  attainment  of  the  project 
goal  are  as  follows: 

1.  Develop  a technique  to  determine 
feasible  fuel  treatment  methods  and  a 
means  to  estimate  their  effectiveness 
and  cost  under  various  sale  conditions. 

2.  Provide  a means  to  examine  the 
influence  of  various  management 
decisions  on  fuel  treatment  costs  and 
effectiveness  (sensitivity  analysis). 

3.  Develop  a scheme  to  quantify  fuel 
treatment  objectives. 

4.  Develop  a computer  program  based 
on  existing  known  relationships  which  can 
easily  be  used  by  a land  manager  to  pre- 
dict fuel  loadings  by  size  classes  prior 
to  harvesting  for  a particular  sale 
(fuel  prediction  model ). 

The  final  results  of  the  project  will  be 
the  development  of  a decision-making 
model  similar  to  the  flow  chart  in 
figure  1.  In  reference  to  figure  1,  a 
decision-maker  is  confronted  with  one  of 
two  situations:  either  select  a fuel 

treatment  system  after  harvesting  or 
select  a fuel  treatment  system  prior  to 
harvesting.  As  mentioned  previously,  we 
are  mainly  concerned  with  the  latter 
situation;  however,  the  final  model  could 
handle  either  situation. 
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SITUATION  NO. 


Figure  l.—Fuel  treatment  selection  model  (a  decisvon-mak%ng  tool) 
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Assuming  situation  No.  2 is  pursued,  the 
land  manager  will  use  the  fuel  prediction 
model  to  estimate  fuel  loadings  resulting 
from  harvesting  based  on  his  particular 
set  of  sale  conditions.  At  this  point 
the  land  manager  has  the  option  of  con- 
ducting a sensitivity  analysis.  He  can 
repeatedly  vary  such  things  as  the  silvi- 
cultural prescription,  falling  of  snags 
or  estimated  damage  of  small  unmerchant- 
able trees  and  observe  their  impact  on 
the  fuel  loading  predicted.  In  other 
words,  he  can  see  how  sensitive  fuel 
loading  is  to  particular  variables. 

Once  a fuel  loading  estimate  is  obtained, 
the  land  manager  must  quantitatively 
specify  his  treatment  objectives  which 
will  be  based  upon  a technique  which  will 
be  developed.  In  addition  to  the  physical 
constraint  of  fuel  loading,  other  vari- 
ables such  as  slope,  residual  spacing  and 
soils  are  considered  in  the  model.  Based 
upon  physical  constraints  and  treatment 
objectives,  feasible  treatment  alterna- 
tives are  determined.  In  other  words, 
the  combination  of  treatment  objectives 
and  physical  constraints  serve  as  a filter 
to  eliminate  obviously  infeasible  treat- 
ments. For  example:  if  the  treatment 

objective  was  to  reduce  fuel  loading  in 
fuels  less  than  4 inches  in  diameter  to 
10  tons  or  less  per  acre  and  present 
loading  in  this  size  class  was  greater 
than  this,  obviously  any  method  of  crush- 
ing would  not  meet  the  treatment 
objecti ve. 

There  are  at  least  three  major  categories 
of  treatments.  First,  there  is  a cate- 
gory that  includes  fuel  treatments  which 
are  conducted  after  harvesting  has  been 
completed.  These  treatments  can  take 
place  either  at  selected  sites  such  as 
landings  or  treat  the  slash  in  place  or 
a combination.  Another  category  con- 
siders treating  the  slash  concurrently 
with  harvesting  by  modifying  the  logging 
system.  For  example:  whole-tree  skid- 

ding and  subsequent  chipping  or  inciner- 
ator burning  tops  at  or  near  landings 
during  the  harvesting  operation  or  using 
a fel 1 er-buncher  to  concentrate  tops 
would  be  considered  a treatment  method. 


Assuming  there  are  additional  costs 
associated  with  this  type  of  modified 
logging  as  compared  to  standard  logging 
methods,  these  additional  costs  would  be 
treated  as  slash  treatment  cost.  The 
final  category  of  fuel  treatments 
includes  subsidies  which  are  provided  to 
entice  the  purchaser  to  remove  certain 
components  of  the  total  fuel  loading. 

For  example:  it  may  be  feasible  to  pay 

a certain  amount  for  the  removal  of  cull 
logs  disregarding  whether  the  purchaser 
makes  a profit  by  selling  the  cull  logs. 
If  the  removal  of  the  cull  logs  signifi- 
cantly reduces  the  total  fuel  treatment 
task,  the  subsidy  used  to  obtain  the 
removal  would  be  treated  as  a slash 
treatment  cost. 

After  feasible  treatment  methods  are 
determined  by  considering  fuel  loadings, 
treatment  objectives  and  physical  con- 
straints, efforts  are  directed  toward 
determining  the  respective  cost  of  each 
treatment.  One  can  use  production  rates 
and  operating  costs  to  estimate  the  total 
cost  for  each  feasible  treatment  method. 

A significant  portion  of  the  project  will 
involve  the  development  of  production 
rates  for  various  fuel  treatments  under 
various  conditions.  The  cost  for  a par- 
ticular treatment  will  then  be  the 
product  of  the  production  rate  times  the 
operating  cost.  It  appears  reasonable 
to  assume  that  once  production  rates  are 
developed,  they  will  remain  valid  unless 
new  technology  is  developed.  Therefore, 
treatment  costs  can  easily  be  updated  by 
revising  operating  costs.  Finally,  once 
treatment  costs  are  obtained  for  feasible 
treatments,  the  treatment  selection  could 
be  based  on  a minimum  cost  since  all 
treatments  which  have  been  considered 
meet  the  treatment  objectives.  However, 
the  decision-maker  has  the  option  of 
rejecting  all  treatments  and  utilizing  a 
sensitivity  analysis  to  obtain  a new  set 
of  feasible  treatments  and  associated 
costs.  For  example,  he  could  change  one 
or  all  of  the  following:  silvicultural 

prescription,  logging  method  or  treatment 
objective.  If  this  route  were  pursued, 
then  a new  set  of  feasible  treatments 
and  costs  would  be  obtained. 
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MAJOR  TASKS  AND  RESULTS  TO  DATE 

The  development  of  the  decision-making 
model  for  the  selection  of  fuel  treatment 
alternatives  is  dependent  upon  the  com- 
pletion of  three  major  tasks  as  illus- 
trated by  figure  2: 

1.  Development  of  fuel  treatment 
production  rates,  limitations,  and 
operating  costs. 


2.  Development  of  a scheme  to  quan- 
tify treatment  objectives. 

3.  Development  of  fuel  prediction 
model . 

The  relative  size  of  each  task  is  illus- 
trated in  figure  2 by  the  size  of  the 
circles,  and  the  approximate  amount  of 
each  task  that  has  been  completed  is 
shown  by  shading. 


Figure  2. — Major  tasks. 
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Fuel  Prediction  Model 

As  shown  by  figure  2,  the  fuel  prediction 
model  has  been  completed.  The  model  con- 
sists of  a user-oriented  computer  program 
that  uses  existing  knowledge  and  data 
that  is  commonly  available  on  Forest 
Service  sales  in  Region  1.  The  model 
utilizes  sale  cruise  data,  Stage  II 
inventory  data  and  certain  information 
provided  by  the  user. 

The  linear  regressions  developed  by 
Fahnestock  for  Region  1 tree  species  are 
used  to  estimate  crown  weights.!/  These 
equations  take  into  consideration  species, 
d.b.h.,  and  crown  length  to  estimate 
crown  weight.  In  addition,  equations 
developed  by  Turnbull  and  Hoyer  to  cal- 
culate tree  volumes  by  utilizing  tariff 
numbers  are  used.!/  Briefly,  the  tariff 
number  system  is  another  way  of  express- 
ing tree  form.  The  tariff  number  is 
simply  the  cubic-foot  volume  of  a tree 
with  I square  foot  of  basal  area.  So 
the  higher  the  tariff  number,  the  better 
the  form  of  the  tree.  The  tariff  number 
system  of  computing  volume  was  used 
because  the  authors  have  developed  equa- 
tions to  calculate  cubic-foot  volumes  in 
terms  of  total  tree,  excluding  the  stump, 
and  also  volumes  to  a 4-,  6-,  and  8-inch 
top.  Therefore,  one  can  readily  estimate 
the  volume  in  the  unutilized  top  stem 
based  on  the  merchantability  standard  of 
the  sale  by  taking  the  difference  between 
the  total  tree  volume  and  the  volume  to  a 
specified  top  size. 

Four  major  sources  of  slash  are  con- 
sidered in  the  fuel  prediction  model: 

1.  Harvested  trees 

a.  Crowns 

b.  Unutilized  top  stems 

c.  Cull  and  breakage 

2.  Snags  that  are  cut  and  not 
utilized. 


3/  Fahnestock,  George  R.  Logging  slash 
flanmability.  Int.  For.  and  Range  Exp.  Sta.  , 
Res.  Paper  No.  58,  pp.  26-33,  May  1960. 

4/  Turnbull,  K.  J.  and  G.  E.  Hoyer. 
Construction  and  analysis  of  comprehensive  tree- 
volume  tariff  tables,  Resource  Mgt.  Rpt.  No.  8, 
State  of  Washington,  Dept,  of  Natural  Resources, 
p.  56,  Aug.  1965. 


3.  Trees  unmerchantable  because  of 
size  that  are  damaged  by  logging. 

a.  Crowns 

b.  Boles 

4.  Existing  fuels  prior  to 
harvesting. 

Figure  3 is  a flow  chart  which  illus- 
trates the  logical  sequences  that  are 
involved  in  the  fuel  prediction  model. 

The  program  is  now  in  the  Region  1 
computer  library  and  is  available  to 
field  units.  Appendix  A is  a set  of 
instructions,  which  has  been  distributed 
to  all  Region  1 field  units,  that 
explains  the  use  of  the  fuel  prediction 
model.  Several  field  units  have  used 
the  program  and  others  plan  to  use  it  in 
the  near  future.  The  purpose  of  the 
fuel  prediction  model  is  to  provide  a 
"ballpark"  estimate  of  resulting  fuels 
prior  to  actual  harvesting.  The  accuracy 
of  the  predictions  is  not  known  at  this 
time;  however,  the  model  incorporates 
the  best  information  that  is  presently 
available  for  fuel  prediction.  As  better 
information  becomes  available  through 
research,  the  model  can  easily  be 
updated  and  improved.  The  reasonable- 
ness of  the  fuel  predictions  of  the  model 
will  be  determined  by  the  experience  of 
field  personnel  using  the  program. 

In  the  instructions  (appendix  A)  is  a 
computer  printout  of  a "run"  using  cruise 
data  and  Stage  II  inventory  data  from  the 
Bear  Run  Timber  Sale  on  the  Missoula 
Ranger  District,  Lolo  National  Forest. 

In  this  particular  case,  the  predicted 
fuel  loading  was  within  10  percent  of 
what  was  actually  measured  on  the  ground 
after  harvesting  had  been  completed. 

The  planar  intersect  method  as  described 
by  Brown  was  used  to  measure  the  fuel 
loading.!/  As  one  can  see,  the  program 
requires  a significant  amount  of  input 
from  the  user  but  in  doing  so  fs  very 
flexible.  A sensitivity  analysis  can  be 
achieved  by  varying  such  things  as: 


5/  Brown,  James  K.  A planar  intersect 
method  for  sampling  fuel  volume  and  surface  area. 
Forest  Science  17(1):  96-102,  1971. 
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Figure 


3. — Computational  diagram  of  fuel  prediction  model. 
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percentage  estimates  of  cull  and  break- 
age of  harvested  trees;  number,  size, 
and  species  of  trees  harvested;  percent- 
age estimate  of  trees  unmerchantable 
because  of  size  damaged  by  logging  and 
so  on. 

Quantify  Treatment  Objectives 

As  depicted  by  figure  2,  no  work  has  been 
completed  in  devising  a scheme  to  quan- 
tify fuel  treatment  objectives.  However, 
it  is  expected  that  slash  treatment 
objectives,  whether  they  are  for  fire 
hazard  reduction,  esthetics,  site  prepa- 
ration, or  whatever,  can  be  expressed 
quantitatively  in  terms  of  fuel  loading 
by  size  classes.  To  estimate  these  fuel 
loadings,  we  plan  to  interview  experts 
in  fire  control,  insect  and  disease  con- 
trol, wildlife  management,  reforestation, 
and  landscape  management.  During  the 
interviews,  photographs  will  be  used  to 
illustrate  various  fuel  loadings.  Based 
upon  these  interviews,  a scheme  will  be 
devised  to  quantify  various  slash  treat- 
ment objectives.  The  scheme  may  appear 
as  follows: 

Treatment  Objective 

1.  Esthetics  - reduce  slash  to 

x tons/ac.  with  nondisturbance 
of  duff. 

2.  Fire  hazard  abatement  - reduce 
to  £ x tons/ac.  plus  reduce  £ 
4-inch  size  class  to  £ y tons/ 
ac.  and  an  average  height  of 
slash  £ z feet. 

3.  Site  preparation  - reduce  slash 
to  £ x tons/ac.  with  y percent 
mineral  soil  exposed. 

4.  Wildlife  - reduce  to  x tons/ac. 
and  remove  all  barriers. 

5.  Other  (e.g..  Insect  and  Disease 
Control)  - reduce  to  £ x tons/ac. 

It  may  turn  out  that  the  final  model  will 
present  the  land  manager  with  a range  of 
values  for  each  objective  and  require  him 


to  select  the  value  which  is  appropriate 
for  his  particular  situation.  Whatever 
technique  is  devised  to  quantify  treat- 
ment objectives,  it  will  not  be  perfect. 
However,  as  land  managers  become  more 
familiar  with  fuel  loadings,  the  tech- 
nique can  easily  be  revised.  The  impor- 
tant point  is  that  the  decision-maker  is 
required  to  think  of  his  objectives  in  a 
quantitative  sense  and  an  attempt  will 
be  made  to  help  him  accomplish  this. 
Refinement  of  a technique  to  quantify 
fuel  treatment  objectives  would  appear 
to  be  a worthwhile  future  research  proj- 
ect by  one  of  the  Forest  Service  experi- 
ment stations. 

Production  Rates,  Limitations,  and 
Operating  Costs 

The  first  step  in  this  part  of  the  proj- 
ect involved  identifying  the  treatment 
methods  to  be  considered.  Appendix  B 
illustrates  those  treatment  methods 
which  will  probably  be  considered.  In 
regard  to  production  rates,  we  hope  to 
develop  relationships  similar  to  those 
illustrated  by  figure  4 for  all  treat- 
ment methods.  Figure  4 is  a hypothetical 
situation  showing  a linear  relationship, 
based  on  a regression  analysis,  between 
production  rates  in  acres/day  and  slope 
for  various  fuel  loadings.  Refer  to 
appendix  C for  a complete  discussion  of 
the  statistical  analysis  which  was  done 
with  hypothetical  data  to  develop  the 
prediction  equation  shown  in  figure  4. 

It  is  recognized  that  there  are  a large 
number  of  variables  which  may  influence 
production  rates.  However,  we  plan  to 
consider  only  those  variables  which 
appear  to  be  most  significant  and  can  be 
readily  identified  such  as  fuel  loading 
and  slope.  Other  variables  such  as 
residual  spacing,  operator  efficiency, 
number  and  height  of  stumps,  soils,  etc., 
will  be  accounted  for  by  making  assump- 
tions in  the  production  schedules  that 
are  developed.  The  production  rates  can 
then  be  adjusted  up  or  down  when  sale 
conditions  are  significantly  different 
from  the  conditions  assumed.  Guidelines 
for  making  these  adjustments  will  be 
developed. 
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Figure  4. — Hypothetical  production  rates. 


Production  rate  data  and  limitations  of 
various  slash  treatment  equipment  have 
been  obtained  from  manufacturers  and 
past  studies  by  the  Forest  Service  and 
others.  In  general,  the  production  rate 
data  are  not  very  informative  primarily 
due  to  variations  in  working  conditions, 
operator  efficiency,  and  accounting 
practices.  Much  of  the  data  will  not  be 
very  useful  in  developing  a slash  treat- 
ment model.  Therefore,  the  bulk  of  the 


data  needed  for  the  development  of  pro- 
duction rates  will  be  generated  through 
the  use  of  expert  opinion.  It  would  be 
impossible  to  generate  the  necessary  data 
through  controlled  time  and  motion  studies 
within  the  time  frame  of  the  project. 
Furthermore,  it  is  doubtful  whether  all 
the  variables  influencing  production  rates 
could  be  integrated  with  time  and  motion 
studies.  Therefore,  it  appears  that  the 
use  of  expert  opinion  is  logical. 
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Expert  opinion  will  be  sought  from  not 
only  the  Forest  Service  but  also 
companies  and  contractors  in  the  region. 
Criteria  will  be  developed  for  selecting 
experts  for  each  treatment  method.  The 
following  factors  will  be  considered: 

1.  Total  experience  with  a treatment 
method  (acres  treated,  number  of  jobs  and 
years). 

2.  Also,  experience  under  various 
sale  conditions. 

3.  Willingness  to  cooperate; 
sincerity. 

Since  most  field  personnel  are  not 
familiar  with  fuel  loadings,  a series 
of  8 x 10  color  photographs  from  the 
National  Fuel  Inventory  and  Classifica- 
tion System  illustrating  various  sale 
conditions  and  fuel  loadings  will  be 
utilized  in  an  attempt  to  supplement 
this  lack  of  experience.  Basically,  each 
expert  will  be  asked  to  estimate  the  pro- 
duction rate  for  a particular  treatment 
method  under  various  sale  conditions.  It 
is  anticipated  that  the  production  rate 
estimates  obtained  from  various  experts 
for  a particular  fuel  treatment  method 
and  a particular  fuel  loading  class  and 
slope  will  not  be  the  same.  Therefore, 
it  may  be  necessary  in  some  cases  to 
attempt  to  reduce  this  variation.  The 
Delphi  Method  will  be  utilized  to  achieve 
the  reduction  in  variation.  Briefly,  the 
Delphi  Method  consists  of  initially 
soliciting  estimates  from  a group  of 
experts  on  an  individual  basis  while 
maintaining  anonymity.  Then  certain 
information  such  as  the  median,  range  of 
estimates  and  the  individual's  estimate 
is  fed  back  to  the  participant  and  he  is 
asked  if  he  desires  to  revise  his  previ- 
ous estimate.  Each  round  or  iteration 
tends  to  draw  the  estimates  closer 
together. 

In  addition  to  seeking  an  estimate  of 
production  rates  for  various  conditions, 
each  expert  will  be  asked  to  estimate  the 
effect  other  variables  such  as  rocks, 
soil,  residual  spacing,  etc.  would  have 
on  his  original  estimate  in  terms  of  a 
percentage  increase  or  decrease.  Finally, 
once  data  has  been  obtained,  we  plan  to 
use  regression  and  correlation  analysis 
to  explore  the  relationship  of  production 


rates  and  various  sale  conditions.  Hope- 
fully, relationships  similar  to  those 
illustrated  by  figure  4 can  be  determined. 

In  February  1974  the  Colville  Ranger 
District,  Colville  National  Forest,  was 
visited  to  discuss  the  feasibility  of 
utilizing  expert  opinion  to  develop  pro- 
duction rates.  Cliff  Lehman,  forestry 
technician  on  the  District,  has  been 
working  with  tomahawks  to  treat  slash 
since  the  early  1 960 1 s and  has  had  con- 
siderable experience  with  dozer  piling. 

To  get  a better  feel  for  some  of  the 
possible  problems  that  one  may  encounter 
in  soliciting  production  rate  estimates 
from  experts,  we  used  Mr.  Lehman  to  make 
a test  run  for  dozer  piling  and  tomahawk- 
ing. Two  variables  were  considered 
directly:  slope  and  fuel  loading. 

Photographs  were  used  to  illustrate 
various  fuel  loadings.  Other  important 
variables  were  considered  indirectly  by 
stating  some  basic  assumptions  and  then 
revising  the  estimate  if  these  assump- 
tions were  not  met.  The  basic  assump- 
tions for  both  dozer  piling  and  tomahawk- 
ing were  as  follows: 

1.  Assume  the  soil  to  be  firm  and 

dry. 

2.  Assume  the  spacing  of  residuals 
to  be  approximately  20  feet  by  20  feet. 

3.  Assume  brush  to  be  approximately 
equivalent  in  density  to  a Douglas-fir/ 
ninebark  habitat  type. 

4.  Assume  stump  heights  are  12 
inches  or  less. 

5.  Assume  there  are  approximately 
25  stumps  per  acre. 

6.  Assume  average  efficiency  of  a 
contractor  operator. 

7.  Assume  an  8-hour  production  day. 

In  reference  to  dozer  piling,  two  more 
assumptions  were  made: 

1.  Assume  that  approximately  75 
percent  of  the  fuel  will  be  piled, 
leaving  25  percent  unpiled. 

2.  Assume  that  approximately  25 
percent  of  the  mineral  soil  will  be 
exposed  (scarified). 
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The  production  rate  estimates  for  vari- 
ous slopes  and  fuel  loadings  provided 
by  Mr.  Lehman  are  illustrated  by 
figure  5.  In  general,  if  one  or  more 
of  the  basic  assumptions  are  not  met 
for  a specific  sale  area,  the  production 
rates  shown  in  figure  5 should  be 
adjusted  accordingly. 

For  slopes  up  to  29  percent: 

1.  If  conditions  are  worse  than  the 
basic  assumptions,  decrease  estimate  up 
to  30  percent  for  the  < 20  tons/acre 
fuel  loading  class  and  up  to  25  percent 
for  other  fuel  loading  classes  depending 


upon  the  severity  of  the  deviation  from 
the  basic  assumption. 

2.  If  conditions  are  more  favorable 
for  production  than  the  basic  assumptions, 
increase  the  production  estimate  up  to  10 
percent  for  all  fuel  loading  classes. 

For  slopes  30  percent  and  greater: 

1.  For  poor  conditions,  decrease  up 
to  20  percent. 

2.  For  favorable  conditions, 
increase  up  to  28  percent. 


Figure  5. — Dozer  piling  (medium-size  D-6). 


In  addition,  the  estimated  production 
rates  should  be  adjusted  as  follows: 

1.  Reduce  production  rates  approx- 
imately 25  percent  if  the  majority  of 
the  fuel  consists  of  full-length  small 
trees  that  were  damaged  by  logging.  This 
type  of  fuel  is  more  difficult  to  pile 
than  normal  logging  slash. 

2.  Reduce  production  rates  approx- 
imately 5 percent  for  each  road  segment 
on  slopes  of  20  percent  or  greater  in 
the  area  to  be  treated. 

In  reference  to  tomahawking,  the  same 
seven  basic  assumptions  were  made.  Pro- 
duction rate  estimates  for  a 3-foot, 
single  tomahawk  mounted  on  a J.D.-450 


or  equivalent  are  illustrated  by  figure 
6.  The  estimates  should  be  decreased 
for  conditions  worse  than  assumed  by  up 
to  30  percent.  Also,  the  presence  of 
rocks  or  rock  barriers  would  cause  a 
decrease  in  production  rates.  Stump 
heights  greater  than  12  inches  would 
also  cause  a decrease.  If  the  slash  is 
arranged  in  concentrations,  a decrease 
should  be  made  since  the  slash  will  need 
to  be  spread  before  treating.  The  pro- 
duction rate  estimates  should  be  increased 
up  to  10  percent  for  conditions  more 
favorable  than  those  assumed.  A 6-inch 
stump  height  is  most  desirable  for  toma- 
hawking, and  the  Colville  Ranger  District 
has  begun  to  require  this  in  their  timber 
sale  contracts  for  areas  that  will  sub- 
sequently be  treated  with  the  tomahawk. 


Figure  6. — Tomahcork  (3- foot  single). 
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In  reference  to  the  limitations  of  dozer 
piling  and  tomahawking  in  partial  cuts, 
Mr.  Lehman  stated  that  the  40  to  60 
tons/acre  fuel  loading  class  was  a maxi- 
mum. Fuel  loadings  beyond  this  range 
should  be  considered  infeasible  since 
the  production  rates  would  be  so  low. 
Maximum  slopes  for  dozer  piling  are 
approximately  35  percent  and  are  approx- 
imately 40  percent  for  tomahawking  and 
would  involve  an  uphill  deadhead  trip. 
Maximum  slope  which  would  facilitate 
both  up  and  down  passes  for  tomahawking 
is  25  percent.  Finally,  tomahawking  is 
not  very  feasible  on  material  larger 
than  4 inches  in  diameter. 

One  problem  anticipated  in  soliciting 
production  rate  estimates  from  experts 
in  regard  to  various  slopes  and  fuel 
loadings  is  the  general  lack  of  knowledge 
or  familiarity  with  fuel  loadings.  Most 
experts  are  accustomed  to  thinking  of 
fuel  in  terms  of  rate  of  spread  and 
resistance  to  control  which  means  differ- 
ent things  to  different  people.  Fuel 
loading  in  terms  of  tons/acre  by  various 
fuel  size  classes  is  a more  precise 
means  of  describing  fuels,  but  few 
people  have  had  much  experience  in  this 
area.  Therefore,  to  bridge  this  general 
lack  of  knowledge,  photographs  of  vari- 
ous fuel  loadings  were  used  when  pro- 
duction rate  estimates  were  solicited 
from  Mr.  Lehman.  Even  though  the  photo- 
graphs were  not  of  the  size  and  quality 
we  plan  to  use  when  actually  interviewing 
experts,  it  is  surprising  how  well  the 
photographs  appeared  to  bridge  the  gap. 
After  viewing  a photograph  illustrating  a 
particular  fuel  loading,  Mr.  Lehman  would 
imnediately  recall  one  or  more  timber 
sales  which  looked  similar  to  the  photo 
before  receiving  treatment.  This 
increased  our  confidence  in  the  use  of 
photos  to  supplement  the  general  lack  of 
knowledge  of  fuel  loadings.  Based  upon 
this  trial  run  with  Mr.  Lehman,  it 
appears  that  obtaining  production  rate 
estimates  by  utilizing  expert  opinion 
is  feasible. 

Operating  costs  regarding  the  use  of  a 
particular  piece  of  equipment  can  be 
developed  by  considering  the  cost  of 
leasing  or  the  cost  of  owning.  The  cost 


of  leasing  when  available  is  a much 
simpler  cost  to  obtain  and  will  be  used. 
It  should  accurately  reflect  the  true 
cost  of  using  the  equipment.  Finally, 
standard  labor  costs  including  fringe 
benefits,  average  overhead  costs  and 
normal  profits  will  be  used  in  develop- 
ing operating  costs  for  a particular 
treatment  method.  The  final  decision- 
making model  will  be  constructed  so 
operating  costs  can  be  adjusted  to  meet 
local  conditions. 


CONCLUSIONS 

1.  Our  problem  analysis  indicated 
that  most  land  managers  have  no  reliable 
or  accurate  method  to  predict  fuel  resi- 
dues resulting  from  a particular  partial 
cut  sale  prior  to  harvesting.  No  method 
exists  for  quantitatively  specifying  fuel 
treatment  objectives  to  determine  feasible 
fuel  treatment  methods.  And,  techniques 
are  not  available  to  estimate  the  associ- 
ated costs  of  a particular  fuel  treatment 
for  a particular  set  of  sale  conditions. 
Therefore,  the  goal  of  the  project  is  to 
develop  a decision-making  model  to  pre- 
dict fuel  loadings,  fuel  treatment  needs, 
treatment  alternatives  and  associated 
costs  prior  to  harvesting  for  a specific 
set  of  sale  conditions.  The  exact  format 
of  the  model  has  not  been  determined  but 
may  consist  of  a computer  program  or  an 
algorithm  utilizing  charts  and  a 

cal cul ator. 

2.  The  fuel  prediction  model  which 
has  been  completed  appears  to  provide  a 
reasonable  method  to  visualize  resulting 
fuels  prior  to  harvesting.  Accuracy  of 
the  model  will  be  determined  by  the 
experience  of  field  personnel  using  the 
tool . 

3.  Accurate  data  for  developing 
fuel  treatment  production  rate  schedules 
which  consider  significant  variables  are 
not  available  even  for  the  most  common 
treatments  and  must  be  obtained. 

4.  A test  on  the  Colville  Ranger 
District  indicated  that  the  use  of 
expert  opinion  to  generate  production 
rate  data  is  feasible. 
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RECOMMENDATIONS 


1.  Continue  the  project  through  fiscal  year  1975. 

2.  Develop  fuel  treatment  production  rate  schedules  and  a scheme 
to  quantify  fuel  treatment  objectives  to  be  used  in  the  final  decision- 
making model  by  utilizing  expert  opinion. 

3.  Monitor  field  use  of  fuel  prediction  model  and  provide  tech- 
nical assistance. 

4.  Forest  Service  research  should  consider  a project  to  develop 
techniques  for  quantifying  fuel  treatment  objectives. 
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United  States  Department  of  Agriculture 

FOREST  SERVICE 

Region  1,  Missoula,  Montana  59801 


APPENDIX  A 


April  19,  1974 


REPLY  TO:  5150 


SUBJECT  Interim  Instructions  for  Use  of  Fuel  Prediction  Model 


Attached  are  Interim  Instructions  for  Use  of  Fuel  Prediction 
Model.  This  program  was  developed  by  Don  Weatherhead  at 
MEDC. 

This  is  a program  that  can  be  used  from  a forest  computer 
terminal  and  should  help  materially  in  making  better  slash 
predictions.  As  stated  in  the  introduction  to  the  program, 
this  model  uses  the  best  information  available  and  as  better 
or  more  information  is  developed  it  will  be  updated. 

The  model  isn’t  difficult  to  use,  but  if  problems  arise  in 
your  program,  instructions  on  who  to  call  are  already  in  the 
computer.  This  program  has  had  limited  testing  and  accuracy 
is  unknown,  but  we  hope  you  will  give  it  a try  and  report 
back  your  results. 

This  is  a long  needed  tool  and  we  have  hopes  that  this,  along 
with  other  forthcoming  information  will  help  you  get  on  top 
of  the  slash  job.  Copies  are  punched  for  insertion  in  the 


TO  Forest  Supervisors 


Division  of  Fire  Management 


cc:  Districts 


Attachment 
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INTERIM  INSTRUCTIONS  FOR  USE  OF  FUEL  PREDICTION  MODEL 


By  DONALD  J.  WEATHERHEAD 
April  18,  1974 


The  Missoula  Equipment  Development  Center  (MEDC)  is  currently  work- 
ing on  a project  involving  fuels  management  on  partial  cuts.  The 
name  of  the  project  is  ED&T  2423  - Fuel  Treatment  Systems  for 
Partially  Cut  Stands.  The  overall  objective  of  the  project  is 
to  develop  a decision-making  tool  to  be  used  by  men  in  the  field 
to  select  desirable  fuel  treatment  systems  in  partial  cuts  prior 
to  harvesting  based  on  the  specific  conditions  of  the  sale.  A 
Project  Record  will  be  available  after  July  1974  which  will  explain 
the  project  in  detail.  One  aspect  of  the  decision-making  tool 
depends  on  the  ability  of  field  personnel  to  predict  fuel  loading 
resulting  from  harvesting  in  terms  of  tons  per  acre  by  size  classes. 
Therefore,  MEDC  explored  the  feasibility  of  developing  a fuel  pre- 
diction model  based  on  presently  known  relationships  and  utilizing 
readily  available  information.  As  a result,  a fuel  prediction 
model  relative  to  Region  1 has  been  developed  in  the  form  of  a 
user-oriented  computer  program.  The  computer  program  is  designed 
to  be  operated  by  a user  at  a remote  computer  terminal.  The  pro- 
gram requests  certain  information  which  is  supplied  by  the  user 
and  then  provides  a predicted  fuel  loading  in  terms  of  tons  per 
acre  by  size  classes.  The  relative  accuracy  of  the  predictions  is 
unknown;  however,  the  program  is  based  upon  the  best  knowledge 
available  regarding  fuel  prediction.  As  better  knowledge  becomes 
available  through  research,  the  model  can  easily  be  updated.  The 
Center  is  very  much  interested  in  the  experiences  field  personnel 
have  in  the  use  of  the  model  regarding  the  relative  accuracy  of 
predictions  and  other  comnents  they  may  have. 

The  fuel  prediction  model  considers  slash  creation  from  six 
sources : 

1.  Crowns  of  harvested  trees. 

2.  Unutilized  portion  of  the  top  stem  of  harvested  trees. 

3.  Estimated  defect  and  breakage  that  will  be  left  as  slash. 

4.  Snags  that  are  cut  and  not  utilized. 

5.  Small  unmerchantable  trees  damaged  by  logging  which  are 
left  as  slash.  This  includes: 
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a.  Crowns 


b.  Boles 

6.  Existing  fuels  prior  to  harvesting  if  known. 

Crown  weights  of  harvested  trees  and  small  trees  damaged  by  logging 
are  computed-|by  using  Fahnestock's  equations  for  tree  species  found 
in  Region  1.—  Tarif  volume  equations  are  used  to  calculate  cubic- 
foot  volumes  of  top  stems,  defect  and  breakage,  snags  and  boles  of 
small  unmerchantable  trees  damaged  by  logging.  The  cubic-foot  vol- 
umes are  then  converted  to  weight  by  considering  the  relative 
specific  gravity  of  each  species. 

Information  needed  to  use  the  fuel  prediction  model  can  be  obtained 
from  the  sale  cruise  and  Stage  II  inventory.  The  model  has  been 
designed  to  be  used  for  estimating  fuels  from  partial  cutting  but 
can  also  be  used  for  clearcuts.  It  is  reconmended  that  fuel  pre- 
dictions be  made  by  timber  types  or  cutting  units  in  partial  cuts. 
The  information  needed  and  its  source  is  as  follows: 

1.  Number  of  acres  in  timber  type 

2.  Including  snags  as  one  species,  the  number 
of  species  to  be  harvested. 

3.  For  harvested  trees  by  species:  Sale  Cruise 

a.  Average  DBH 

b.  Average  height 

c.  Total  number  of  trees  harvested  in 
timber  type. 

4.  For  trees  unmerchantable  because  of  size: 

a.  Number  of  diameter  classes  in  timber  type 
(Consider  all  trees  in  a diameter  class  as 
one  species.) 

Stage  II 

b.  By  species: 

(1)  Average  DBH 

(2)  Number  of  trees  per  year 


V Fahnestock,  George  R.,  Logging  Slash  FI amnabi 1 i ty . Inter- 
mountain Forest  and  Range  Experiment  Station,  Research  Paper  No.  58, 
pp.  26-33,  May  1960. 
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5.  Other  information  needed: 


a.  Estimate  of  defect  and  breakage 
left  as  slash. 


Based  on  user's 
estimate 


b.  Estimate  of  percentage  of  trees 
unmerchantable  because  of  size 
that  are  expected  to  be  damaged 
by  logging  and  left  as  slash. 


c.  Size  of  minimum  merchantable  top 
of  harvested  trees. 


Timber  Sale 
Contract 


d.  Fuel  loading  prior  to  harvest- 
ing. 


Preharvesting  Fuel  In- 
ventory if  available. 


e.  Tarif  number  by  species 


Determined  by  user 


Tarif  numbers  have  been?raentioned  in  the  preceding  discussion  with 
no  explanation  of  them.-  Briefly,  a tarif  number  is  simply  the 
cubic-foot  volume  of  a tree  with  one  square  foot  of  basal  area. 

For  example,  a tarif  number  of  27  means  that  a tree  with  one  square 
foot  of  basal  area  (approximately  13.5"  DBH)  will  have  a volume  of 
27  cubic-feet.  Therefore,  the  larger  the  tarif  number,  the  more 
volume  found  in  a tree.  So  the  tarif  number  is  really  another  way 
of  expressing  the  form  of  a tree.  Tarif  numbers  are  used  in  the 
fuel  prediction  model  because  cubic-foot  volume  can  easily  be 
computed  using  equations  for  the  entire  tree  including  the  top 
and  volume  to  a 4-,  6-,  and  8-inch  top.  Therefore,  it  is  simple 
to  estimate  the  volume  in  the  unutilized  portion  of  the  top  stem 
in  harvested  trees.  For  example,  if  the  minimum  top  DIB  is  6 
inches,  one  can  determine  the  volume  in  the  entire  tree  including 
the  top  and  then  determine  the  volume  to  a 6-inch  top.  The 
difference  between  the  two  volume  calculations  is  the  volume  in 
the  unutilized  top  stem. 

Tarif  numbers  can  be  determined  by  using  a computer  program  called 
*TARIF  which  is  stored  in  the  CSC  computer  in  Los  Angeles. 

Region  1 users  have  access  to  this  computer  facility.  Exhibit  1 
is  a photocopy  of  a printout  of  a "run".  As  with  the  fuel  pre- 
diction model,  this  program  also  depends  on  user  supplied 
information.  Tarif  numbers  should  be  determined  for  each  species 
by  timber  type.  As  noted  in  Exhibit  1,  the  program  requests 
the  average  DBH  and  average  height  by  species  and  then  the 


2/  Turnbull,  K.J.  and  G.E.  Hoyer,  Construction  and  Analysis 
of  Comprehensive  Tree-Volume  Tarif  Tables,  Resource  Management 
Report  No.  8,  State  of  Washington,  Dept,  of  Natural  Resources, 
p.  64,  April  1965. 
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tarif  number  is  computed.  The  program  is  designed  to  accept  DBH's 
to  the  nearest  one-inch  class  for  diameters  5-10  inches  and  the 
nearest  even  two-inch  class  for  larger  diameters. 

Once  tarif  numbers  and  other  information  have  been  obtained,  the 
user  can  make  a fuel  prediction.  The  fuel  prediction  model  is 
called  *FUEL  and  is  also  stored  at  the  CSC  computer  center. 

Exhibit  2 is  a photocopy  of  a "run"  using  data  from  the  Bear  Run 
Timber  Sale  on  the  Missoula  Ranger  District.  Referring  to 
Exhibit  2,  the  first  half  or  more  of  the  program  consists  of 
instructions  and  questions  to  obtain  information.  The  last  line 
of  the  program  gives  the  predicted  fuel  loading  by  size  classes. 
Results  of  preceding  computations  of  various  fuel  components 
are  also  printed  out  prior  to  the  final  fuel  loading  estimate. 
These  results  may  be  of  value  to  some  users.  Also,  it  should  be 
noted  that  the  program  is  designed  to  be  very  flexible  and  a 
sensitivity  analysis  can  easily  be  accomplished.  For  example, 
one  can  explore  the  effect  upon  fuel  loading  by  varying  such 
things  as: 

1.  Percentage  estimate  of  trees  unmerchantable  because 
of  size  damaged  by  logging. 

2.  Size  and  number  of  harvested  trees. 

3.  Minimum  top  size. 

4.  Snags  (number  cut  and  not  utilized). 

5.  Estimated  defect  and  breakage  left  as  slash. 

As  mentioned  earlier,  MEDC  is  interested  in  the  experiences  of 
field  personnel  using  the  fuel  prediction  model  and  would  appre- 
ciate their  opinion  regarding  the  accuracy  and  usefulness  of  the 
model.  Finally,  it  should  be  restated  that  the  model  is  an 
attempt  to  utilize  existing  knowledge  and  readily  available 
information  to  predict  fuel  loading.  The  accuracy  of  the  model 
is  unknown.  The  total  fuel  loading  estimate  may  be  reasonably 
accurate  as  well  as  the  4-inch  and  greater  size  class.  However, 
the  breakdown  of  size  classes  less  than  4 inches  is  based  on  a 
minimal  amount  of  data. 
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EXHIBIT  1 


SAMPLE  "RUN"  OF  TARIF  PROGRAM 


*TARIF  1£:  5?  04.-"  1 74 

THIS  PROGRAM  DETERMINES:  THE  RPPROPRIRTE  TFlRIF  NUMBER  FOR  EACH 
TREE  SPECIES  TO  USE  WITH  THE  FUEL  PREDICTION  MODEL.  YOU 
SHOULD  DETERMINE  TfiRIF  NUMBERS  BY  TIMBER  TYPES  FOR  EACH  SPECIE 
1 HE REF ORE?  MAKE  SEPARATE  ’ RUNS’  WITH  THIS  PROGRAM  FOR  EACH 
TIMBER  TYPE. 


SPECIES:  CODE 

WESTERN  WHITE  PINE  1 

WESTERN  REBCEDAR  £ 

WESTERN  HEMLOCK  3 

ENGELMANN  SPRUCE  4 

GRAND  FIR  5 

SUEHLP I ME  F IF'  6 

LODGEPOLE  PINE  7 

P L i I I D E R U S f i P 1 N E 

DOUGLAS-FIF:  0 

WESTERN  LARCH  10 


HOW  MANY  SPECIES  DO  YOU  WANT  TARIF  NUMBERS  FOR  IN  THIS  TIMBER 
TYPE? 

■~iC 


TO  DETERMINE  TARIF  NUMBERS  BY  SPECIES  FOR  THIS  TIMBER  TYPE? 
USE  THE  ABOUE  NUMERICAL  CODES  FOR  SPECIES  AND  TYPE  IN  THE 
FOLLOWING  DATA  (SEPARATE  EACH  DATA  ENTRY  WITH  A COMMA): 
SPECIES?  HUE  DBH  AND  RUE  HEIGHT. 

: 1 S ? Is?  6 y ? 1 y ? c‘ G ? 1 y y ? 7 ? 1 u ? y y ? y > IS?  OS  ? h ? It'?  t> S 

RUE  AUE  TAR 
SP  DBH  NT  NO 
S 1 L SU  c!  i 

1G  £S  100  £8 

7:  10  58  £9 

S 1 6 56  1 9 . 5 

9 1 6 66  £ 1 . 5 

MOW  AT  17'S 

SRU ’ S : 8 . 3 
READY 
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EXHIBIT  2 


SAMPLE  "RUN"  OF  FUEL  PREDICTION  MODEL 


RUN 

*FUEL.  13:@0  84--'l£.-''74 

THIS  PROGRAM  PREDICTS  FUEL  LOADING  RESULTING  FROM  PARTIAL 
CUTTING  FOR  EACH  TIMBER  TYPE  IN  YOUR  SALE.  THE  PREDICT  I ON 
IS  BASED  UPON  CRUISE  AND  STAGE  II  DATA  WHICH  IS  SUPPLIED 

by  the  user. 

WHEN  DATA  IS  REQUESTED?  USE  THE  FOLLOWING  NUMERICAL  CODES  FOR 
TREE  SPECIES: 


SPECIES:  CODE 

WESTERN  WHITE  PINE  1 

WESTERN  REDCEDAF:  £ 

WESTERN  HEMLOCK  3 

ENCELMANN  SPRUCE  4 

GRAND  FIR  5 

SUBALP I ME  FIR  6 

LODGEF'OLE  PINE  7 

F' U NdERU S A F‘  I N E W 

D u U GLfiS— F IR  9 

WESTERN  LARCH  10 

SNAGS  1 1 


INCLUDING  SNAGS  AS  ONE  SPECIES?  HOW  MANY  COMMERCIAL  SPECIES 

ARE  ! HERE  IN  THIS  l I MbEF:  TYF'EV 


ENG  THE  ABOUE  NUMERICAL  CODES?  TYPE  IN  THE  FOLLOWING  DATA 
FOR  EACH  SPECIES  (SEPARATE  EACH  DATA  ENTRY  WITH  A COMMA): 
SPECIES?  RUE  DEH?  RUE  HEIGHT?  DEFECT?  TOTAL  TREES  PER 
SPECIES?  AND  TAR IF  NUMBER.  REMEMBER  TO  ENTER  ZERO’S  WHERE 
APPROPRIATE  AND  INCLUDE  SNAGS.  DEFECT  REFERS  TO  DEFECT 
AND  BREAKAGE  THAT  YOU  ESTIMATE  WILL  BE  LEFT  AS  SLASH  FUEL. 

Yb  ? i E ? bo  ? S . S ? 11  bO  ? E7  ? 1 0 ? Eb  ? 1 00  ? O ? Eb  ? Eo  ? 7 ? 1 0 ? S S ? 4 . b ? Eb47  ? E'? 

?S ? lb? 5b ? 3 . 3 ? 83£ ? 19.5? 9 ? lb? bb ? 9 . 1 ? S170?  £1.5? 11? 19?  73 ? 1 OO ? 41b?  £ 
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VO  UR  DATA  CONSISTS  OF  THE  FOLLOWING: 


HUE  HUE  i CTRL  T RE 

SF'  DEH  HI  i>EF  TREES  Nu 
6 iE  bo  . y 1 1 66  E 

i El  c!fc>  1 06  U ci6  EE 


10 


4.6 


.64  r 


i i 


19. 


66  9 .1  61 70  £1.5 

78  106  416  £8 


IS  THIS  BflTFl  CORRECT  (YES=1*  NO=Gj? 

?1 

GOOD ! 

HOW  MANY  RCRES  FIRE  IN  THIS  TIMBER  TYPE? 

166 

ns  I ND I CRT  ED  ON  THE  STFlGE  II  INUENTORY?  HOW  MANY  DIRMETER 
CLASSES  ARE  REPRESENTED  IN  TREES  UNMERCHANTABLE  BECAUSE 
OF  SIZE? 

MEET  TYPE  IN  THE  FOLLOWING  DATA  FOR  UNMERCHANTABLE  TREES 
BECAUSE  OF  SIZE  : 

SPECIES?  DBH?  TREES--ACRE?  TAR  IF  NO.  USE  THE  PRECIOUSLY 
STATED  CODE  FOR  SPECIES. 

YOUR  DATA  CAN  BE  OBTAINED  FROM  THE  STAGE  II  INUENTORY. 
ASSUME  ALL  TREES  IN  EACH  DIAMETER  CLASS  ARE  ONE  SPECIES. 


YOUR  DATA  FOR  TREES  UNMERCHANTABLE  BECAUSE  OF  SIZE  IS 
AS  FOLLOWS: 

TREES  TAR IF 
SF'  DBH  ••'AC  NU 
cl  cl  U 4 U cl  i 

6 4 60  £? 


IS  IHIS  DAT  hi  CORRECT  I YES  — 1 ? NU— UJ  Y 
?1 
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THERE  ARE  SEUERRL  ASSUMPTIONS  MADE  IN  THIS  PROGRAM  WHICH 
CAN  BE  ALTERED  IF  YOU  DISAGREE.  FIRST  OF  ALL  IT  IS  ASSUMED 
THAT  APPROXIMATELY  50".  OF  THE  SMALL  UNMERCHANTABLE  TREES 
WILL  BE  DAMAGED  BY  LOGGING  AND  WILL  CONTRIBUTE  TO  THE  SLASH 
PROBLEM . 

DOES  THIS  APPEAR  TO  BE  A REASONABLE  ASSUMPTION  (YES=ls 
NO-0j  ? 

?0 

TYPE  IN  YOUR  PERCENTAGE  ESTIMATE. 

■ i 'i  c 

C---1 

ANOTHER  ASSUMPTION  THAT  HAS  BEEN  MADE  IS  THAT  APPROXIMATELY 
10X  OF  THE  TOTAL  LIMBS  ON  ALL  COMMERCIAL  SPECIES  ARE  DEAD. 
DOES  THIS  APPEAR  TO  BE  REASONABLE  ( YES=  1 .-  N0=8  j ? 

?i 

THIS  PROGRAM  CONSIDERS  A MERCHANTABLE  TOP  DIB  OF  4s 6 OR  S 
INCHES. 

TYPE  IN  THE  MERCHANTABLE  TOR  SIZE  FOR  THIS  SALE  THAT  IS 
NEAREST  TO  ONE  OF  THE  THREE  SIZES  ABOUE. 


IS  THERE  A PREHARUESTING  FUEL  INUENTORY  AUAILABLE  FOR  THIS 
TIMBER  TYPE  fYES=is  NO=0j? 

?1 

TYPE  IN  THE  F UEL  LOADING  ( TONS - ACRE. J FOR  THE  UARIOUS  SIZE 
CLASSES  IN  THE  FOLLOWING  ORDER:  0-  1,-4  INCHES.-  1/4  - 1 IN.  - 

1 - 3 IN.-  3+  IN.  AND  TOTAL. 

73 . 1 - 1 . 5 - 3 . 3 - 4 . 3 - 11. 3 


CROWN 

INFORMATION  BY 

SPECIES 

SPECIES 

CROWN  LENGTH 

CROWN  WEIGHT 

TOTAL  WEIGHT 

n 

44 . 94 

144.937 

184936. 

10 

7 1 . 6 

346. 433 

9909 . 4 

- 

30 . 

75. 4736 

3 1 9 r 5 r . 

45 . 6 

o 1 0 . 3 r 4 

333963 . 

S1 

45.6 

305.574 

1 . 39533E+6 

C R U W N W EIGHT  B Y 3 1 ZE  ULA S S E S 
(TOTAL  POUNDS  FOR  TIMBER  TYPE) 

1/4-1  IN  1-4 IN  TOTAL 

6 1 4943 . 7 1 9069 . 3 . 0937SE+6 


O- 1 /4 1 N 

!*  ? r »'  d m 
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UNMERCHANTABLE  TOP  WEIGHTS 
(TOTAL  POUNDS  FOR  TIMBER 
SPECIES  TOR  WEIGHT 

6 48531 . 9 

HI  I 8 1 9 . 1 

r i 8 1 Go  f*  . 


1 


4518. 

1:4954, 


BY  SPECIES 
TYPE) 


WEIGHT  OF  TOP  STEMS 
(TOTAL  POUNDS  FOR  TIMBER  TYPE) 
- 4- INCHES  TOTAL 

11 1730.  367310. 


WEIGHT  OF  SNAGS 
(TOTAL  POUNDS  FOR  TIMBER  TYRE) 
453183. 


WEIGHT  FROM  DEFECT  AND  BREAKAGE 
(TOTAL  POUNDS  FOR  TIMBER  TYPE) 
58697 1 . 


SMALL  UNMERCHANTABLE  TREES  DAMAGED  BY  LOGGING 
CROWN  WEIGHT  BY  SIZE  CLASSES 
POUNDS.- 'ACRE 


O- L-4  IN 
3439.83 


1.--4-1IN  1-4  IN  TOTAL 

1493.18  1759.76  6691.11 


BOLE  WEIGHTS  OF  TREES  DAMAGED  BY  LOGGING 
POUNDS.--ACRE 
-4- INCHES  lOTAL 

8084. 1 8 8884. 1 8 


FUEL  LOADING  FOR  TIMBER  TYPE  BY  SIZE  CLASSES 

T Of  iS.-ACRE 

0-1, -"4  IN  1 .-"4-1  IN  1-4  IN  4+  IN 

6. lUoUo  4. 09888  7. 99435  8. 88808 

NOW  AT  1 606 
SPINS:  1 . 6 
READY 


TOTAL 

LL  *1*  • *"t  L 
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APPENDIX  B 


POSSIBLE  TREATMENTS  TO  BE  CONSIDERED 


A.  Dozer  Piling 

1.  light  ) 

2.  medium  ) based  on  horsepower  and  weight 

3.  heavy  ) 

B.  Hand  Piling 

C.  Chipping 

1.  Inplace 

a.  light  ) 

b.  medium  ) based  on  size  of  chipper  head  and  horsepower 

c.  heavy  ) 

2.  Selected  sites 

a.  light 

b.  medium 

c.  heavy 

d.  include  skidding  or  bunching  costs 

e.  include  cost  of  disposal  of  chips  or  selling  value 
whichever  is  appropriate. 

D.  Tree-length  skidding  to  landing 

E.  Skidding  Tops 

F.  Masticating  Inplace  (includes  Tomahawk,  Morden  Chopper,  Tree-Eater,  etc.) 

G.  Additional  Cost  of  Using  a Feller-Buncher  to  Concentrate  Tops  for 
Skidding  or  Treatment 

H.  Dozer  Trampling 

I.  Pit  or  Incinerator 

J.  Burying 


-25- 


' 


1 


APPENDIX  C 


J H I’ I Y TO: 
St  1 11.1 1 ( I 
TO: 


1 


UnITEDSTATES  DEPARTMENTOF  AGRiCULTURE 
. FOREST  SERVICE  „ 

Equipment  Development  Center 
Fort  Missoula 
Missoula,  Montana  59801 


7120  Equipment  Development  and  Test  June  27,  1974 


ED&T  2423  - Fuel  Treatment  Systems  for  Partially  Cut  Stands 


A.  H.  Jukkala,  Project  Leader 


The  following  is  a statistical  analysis  of  hypothetical  production 
rate  data.  The  objective  of  this  analysis  is  to  become  familiar 
with  the  procedures  I intend  to  use  with  the  real  data  and  also 
identify  any  statistical  problems  that  could  be  avoided  when  obtain- 
ing the  data.  Exhibit  1 illustrates  the  type  of  data  that  is 
expected  to  be  obtained  for  dozer  piling  slash.  The  data  repre- 
sents production  rates  in  terms  of  acres/day  that  can  be  treated 
with  a medium  size  dozer  (D6)  under  various  slope  and  fuel  loading 
combinations.  The  production  rate  estimates  underscored  are  actual 
estimates  from  an  expert  I interviewed  in  February  1974.  The  other 
estimates  in  each  slope  and  fuel  loading  category  represent  hypo- 
thetical estimates  from  six  other  experts. 

Figures  1-3  illustrate  the  data  plotted  for  three  fuel  loading 
classes  and  various  slopes  ranging  from  0 % to  35%.  Since  the 
plan  is  to  fit  a curve  to  this  data  using  regression  analysis,  some 
assumptions  must  be  made.  It  is  assumed  that  all  Y's  (production 
rates)  are  normally  distributed  and  have  the  same  variance  for  all 
values  of  X (slope  and  fuel  loading).  The  plotted  data  on  figures 
1-3  suggests  that  possibly  the  Y's  do  not  have  homogeneous  variance 
since  the  range  in  production  rate  estimates  is  relatively  large 
for  small  values  of  slope  and  fuel  loading  and  decrease  as  slope 
and  fuel  loading  increase.  Therefore,  I attempted  to  compensate 
for  this  possible  heterogeneity  in  variance  with  a logrithmic  trans- 
formation of  the  Y values  as  illustrated  by  figure  4.  However,  as 
the  figure  depicts,  the  possible  heterogeneity  was  reversed  with 
the  transformation . This  suggested  that  possibly  it  was  reasonable 
to  assume  that  the  Y's  did  have  homogeneity.  In  a further  effort 
to  determine  whether  it  was  reasonable  to  assume  homogeneity,  I 
used  Bartlett's  test  as  described  by  Snedecor.-  I wrote  a computer 
program  to  perform  the  calculations  using  Snedecor' s data  for  a test. 


-^Snedecor,  G.  W.  1962.  Statistical  Methods.  Chapter  10, 
pp.  285-287.  Iowa  State  University  Press,  Ames,  Iowa. 


-26- 


2. 


The  program  and  output  are  shown  in  exhibit  2.  I then  used  the  pro- 
gram to  calculate  a Chi  Square  for  my  data  for  each  fuel  loading. 

With  each  fuel  loading  I had  seven  samples  (values  of  slope)  with 
seven  observations  each  (production  rate).  Exhibits  3-5  show  the 
calculated  Chi  Squares  obtained  which  range  from  1.26  to  3.9.  In 
each  case  I was  testing  the  hypothesis  that  the  variance  of  pro- 
duction rate  is  the  same  for  all  values  of  slope.  The  null  hypothesis 
is  as  follows: 


H 


0 • 


2 2 2 2 2 2 2 

Si  = S2  = S3  = S4  = S5  = S6  = S 7 


The  tabulated  value  of  Chi  Square  with  six  degrees  of  freedom  at  the 
95%  confidence  level  is  12.59.  Since  all  of  the  calculated  Chi 
Squares  are  less  than  12.59,  I could  not  reject  the  hypothesis  at 
the  95%  confidence  level.  Therefore,  it  appears  reasonable  to  assume 
that  the  Y's  have  homogeneous  variance. 

I used  a step-wise  multiple  regression  canned  program  from  Computer 
Science  Corporation  in  Los  Angeles  to  determine  a best  fitting 
equation  for  the  data.  The  plotted  data  suggested  that  a curve  may 
fit  better  than  a straight  line.  Therefore,  I considered  the  follow- 
ing independent  variables: 

Xj  = slope  (%) 

x2  = fuel  loading  (tons/ac.) 

X 3 = (Xl) 

X4  = (x2)2 

x5  = (XjX2) 

I began  by  considering  all  of  the  variables  in  the  regression.  The 
model  was  as  follows: 


y = b0  + b1x1  + b2x2  + b3x3  + b4x4  + b5b5 

The  computer  printout  (see  exhibit  7)  suggested  that  regression  due 
to  X2  was  insignificant;  therefore,  I tested  the  hypothesis  that 
b2  = 0 at  the  5%  significance  level  using  an  F-test  in  an  analysis 
of  variance  (exhibit  6).  The  H0:  b2  = 0 could  not  be  rejected  at 
the  5%  level.  This  same  conclusion  was  also  arrived  at  using  the 
calculated  t value  for  b2  from  the  computer  printout.  Therefore, 

I considered  the  optimum  model  to  be: 

y = b0  + biXi  + b3x3  + b4x4  + bsx 


5 5 
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3. 


I then  tested  the  significance  of  each  variable  using  an  F-test 
(exhibit  6).  It  turned  out  that  all  of  the  variables  were  signifi- 
cant at  the  5%  significance  level.  The  final  prediction  equation 
is : 

y = 10.4141  - .298141  (x1 ) + .00190705  (x3)  - .0024816  (x4)  + 
.00231329  (x5)  with  R2  = .753. 

In  other  words,  the  equation  explains  about  75%  of  the  variation  in 
Y.  Or  the  use  of  slope  and  fuel  loading  in  the  equation  accounts 
for  approximately  75%  of  the  variation  in  production  rates. 

Finally,  figure  5 illustrates  the  prediction  equation  for  various 
values  of  slope  and  three  levels  of  fuel  loading.  It  is  interest- 
ing to  note  that  the  interaction  of  slope  and  fuel  loading  is  the 
most  significant  variable.  In  addition  the  relationship  of  produc- 
tion rate  to  slope  and  also  to  fuel  loading  is  curvilinear.  The 
prediction  equation  (curved  line)  was  also  drawn  through  the  plotted 
data  in  exhibits  1-3. 


DONALD  J.  WEATHERHEAD 
Forester 
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EXHIBIT  1 


DOZER  PILING  - MEDIUM,  D6 


Estimates  of  Production  Rates  (acres/day)  by 
Slope  and  Fuel  Loading  Classes 


Tons/acre 

Tons/acre 

Slope 

<20 

20-40 

40-60 

Slope 

<20 

20-40 

40-60 

15 

12 

6 

7 

6 

3 

12 

10 

5 

6 

5 

2.5 

0 

11 

8.5 

4.5 

25 

5.5 

4 

2 

10 

_8_ 

A 

5 

3.5 

1 

9.5 

7.5 

4 

4 

3 

.8 

9 

7 

3.5 

3 

2 

.7 

8 

6 

3 

2.5 

1 .8 

.5 

10 

9 

5 

6 

5 

2.5 

9 

8 

4 

5 

4 

2 

10 

8 

7 

7.5 

2. 

3.5 

3 

30 

4.5 

4 

3.5 

3 

1 

.8 

6.5 

6 

2.8 

3.5 

2.5 

.6 

6 

5 

2.5 

3 

2 

.4 

5 

4 

2 

2.1 

1.5 

.3 

9 

8 

4.5 

5 

4 

2 

8 

7 

3.5 

4 

3 

1.5 

15 

7 

6 

3 

3.5 

2.5 

1 

6 

£ 

2.5 

35 

3 

2 

.5 

5 

4 

2.3 

2.5 

1.8 

.4 

4.5 

3.5 

2.0 

2 

1 .5 

.3 

4 

3 

1 .5 

1 .8 

1 .1 

.2 

8 

6.5 

3.5 

7 

6 

3 

20 

6 

5 

2.5 

5.5 

A 

2 

5 

3 

1 .8 

4 

2.5 

1.5 

3 

2 

1 
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HOMOG 1 11=07 

HTL  TUE  U6--''l3--‘''74 

4 SAMPLES  WITH 

5 OBSERUATIONS  EACH. 

SAMPLE 

SUM  OF  SQUARES  MEAN  SQUARE 

1 

478  1 1 8 

U_ 

396  99 

616  154 

'T 

164  41 

EXHIBIT  2 


LOG  OF  MEfiN  SQUARE 
£. 071 SS 
1 . 99564 
£. 1 875£ 

1.61 £78 


MEfiM  OF  St£  = 103  LOG  OF  MEAN  ST£  = £.01 £84 
DIFFERENCE  = . 1835£7 
UNCURRELTED  CHI  SQUARE  = 1 . 69U36 
DEGREES  OF  FREEDOM  = 3 

CORRECTED  CHI  SQUARE  = 1.53089 


RUNNING  TIME:  .5  SECS  I--0  TIME 


. 4 


SECS 


READY 

LIS7TIH 

1 REM  COMPUTATION  OF  BARTLETT’S  TEST  FOR  HOMOGENEITY  OF  UARIANCE. 

£ REM  SAMPLES  EQUAL  IN  SIZE.  REFER  TO  SNEDECOR  PAGE  £86. 

3 REM  ENTER  DATA  STARTING  WITH  LINE  £Q  - £0O.  FIRST  ENTER  A AND 

5 REM  N.  A=  NUMBER  OF  SAMPLES  ? N=  NUMBER  OF  OBSERUATIONS  PER 

7 REM  SAMPLE.  DATA  IS  ENTERED  FOR  THE  PROBLEM  SNEDECOR  DISCUSSES. 

£0  DATA  4.-5 

30  DATA  40? £4? 46? £Q? 35 
40  DATA  £9? £7? £0? 39? 45 
50  DATA  1 1 ? 3 1 ? 1 7 ? 37  ? 39 
60  DAT  A 1 7 ? 8 1 ? £8  ? 33 ? 8 1 
£10  READ  A?  N 

£15  PRINT  A? "SAMPLES  WITH" ? NS "OBSERUATIONS  EACH." 

£18  PR I NT "SAMPLE"? "SUM  OF  SQUARES" ? "MEAN  SQUARE" ? "LOG  OF  MEAN  SQUARE" 

££0  FOR  J = 1 TO  A 

££5  LET  XI  = Q 

££7  LET  X£  = Q 

£30  FOR  I = 1 TO  N 

835  READ  X 

£40  LET  XI  = XI  + X 

£50  LET  X£  = X£  + XT£ 

£60  NEXT  I 

£6£  LET  X3  = 0 

£64  LET  X4  = O 

£66  LET  X5  = 0 

£70  LET  X3  = X£  - Xli  £.-  N 

£80  LET  X4  = X3---''  C N- 1 i 

£90  LET  X5  = CLG(X4) 

300  LET  X6  = X6  + X4 
310  LET  X7  = X7  + X5 
330  PR  I NT  J ? X3 ? X4 ? X5 
340  NEXT  J 
350  PRINT 
355  LET  XS  = X6--A 
358  LET  X9  = CLGiXS) 

360  PRINT  "MEAN  OF  STS  =" 5 X8"  LOG  OF  MEAN  STS  ="?X9 

370  LET  Y = A*X9  - X7 

380  PRINT  "DIFFERENCE  ="?Y 

390  LET  Y 1 = £.30£6*(N-1)»SY 

400  PRINT  "UNCORRECTED  CHI  SQUARE  ="Y1 

410  PRINT  "DEGREES  OF  FREEDOM  ="?A-1 

4£0  LET  C = 1 + ( CA+1  !i  (3ffiAffi  (N— 1 j ) H 

430  LET  Y£  = Y1--C 

440  PRINT 

450  PRINT  "CORRECTED  CHI  SQUARE  ="?Y£ 

999  END 
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£0 

DATA 

73  7 

30 

DATA 

15? 1 £ ? 11’ 1 0 ’ y . 5 ’ 

9 ’ 8 

40 

DAT  A 

1 0 ’ A ’ 8 1 7 ’ 6 . 5 ’ 6 ’ r 

50 

DAT  A 

9 ! 0 3 7 3 6 3 5 3 4 ■ 5 3 4 

60 

DAT  A 

:~i  t*  cr:  cr  nr  .•« 

m 3 i 3 t.  3 ._i  q ._!  3 _(  3 ‘-j-  3 £ i 

70 

DATA 

7 ? 6 ? 5 . 5 j .■  ’ 4 j 3 ? £ . 

irr 

88 

DATA 

I~i  H !~il  4 . 4 *»  4 „ 4 * - **I  1 

£ . 1 

90 

DATA 

9 3 4 3 4 .03  9 3 «*'  a 4 3 9 3 

1 . 8 

Kur 

\ 

EXHIBIT  3 
10  TONS  / ACRE 


HuMOGl 


11:19  RTL  TUE 


06--'  1 S--74 


SAMPLES  WITH  7 OBSERVATIONS  EACH. 


SAMPLE 

1 


SUM  OF  SQUARES  MEAN  SQUARE 


3£ . 357 1 
18.3571 
£8 . 9£86 
17.5 
1 •!'  a 9 c!  96 
10. 1086 
i u 94957 


5-  99990 
3. 05952 
3.4881 
£.91 667 
£ . 6547 6 
1 . 6847 6 
1 . 308 1 


LOG  OF  MEAM  SQUARE 
.731819 
. 485654 
. 5485 


. 4648 
. 4848 

. 1 1 66 


C.C- 

39 

-:9 


MEAM  OF  STS  = S.9S9S5 
DIFFERENCE  = .£75145 
UNCORRECTED  CHI  SQUARE  - 3.8813 
DEGREES  OF  FREEDOM  = 6 


LOG  OF  MEAM  STS 


CORRECTED  CHI  SQUARE  = 3.57436 


38  DATA 
48  DATA 
58  DATA 
6U  jJHT  A 
78  DAT  A 
30  DATA 
90  DATA 
RUM 


1 £ j 1 0 ? 8 . 5 ? 8 j 7 . U j 7 ’ 6 
y s cl?  7.5;  7 ? 6 ? 5 ? 4 

9 3 •'  3 3 3 5 3 *4'  3 5 ■ 9 3 0 
9 c 5 ? t'  3 9 3 4 ? 3 9 a 5 3 9 3 

9 3 5 3 4 3 9 . 5 3 9 3 9 3 1 o S 

9 3 'T  3 0 . -Ii  ? 9 3 Cl  a ~J  3 9 3 1 . 5 

4 ? 8 ? £ . 5 ? £ ? 1.8?  1 . 5 j i . i 


EXHIBIT  4 
30  TONS/ ACRE 


HOMOG1 


11 


ATL  TUE 


U6--''  1 8.--7  4 


7 SAMFLES  WITH*  7 0ESERUATI0NS  EACH. 


SAMPLE 

SUM  OF  S 

QUARES  MEAN  SQUARE 

LOG  OF  MEAN 

1 

£4. £143 

4.83571 

. 60598 

C. 

18. 3571 

3. 05958 

. 485654 

9 

£8.9886 

3. 4881 

nr  y •“  rr 
° JTil-JOu 

4 

18.3571 

-~j  «~jcr  i~icr~i 

O s i-J-J  Z-  •_»  C. 

. 485654 

ET 

14.8486 

£. 34143 

. 369481 

0 

8 . r 1 4£9 

1 . 45838 

. 168081 

r 

5. 33489 

. 97838 1 

~ 1 . £ 1 636E — £ 

MEAN  OF  STS  = 
DIFFERENCE  = 

£ . 6£986 
. 3003 1 9 

LOG  OF  MEAN  STS  = 

. 4 1 9933 

UNCORRECTED  CHI  SQUARE 
DEGREES  OF  FREEDOM  = 6 

= 4. 14903 

CORRECTED  CHI 

SQUARE  = : 

:L  90138 

35 


' 


EXHIBIT  5 
50  TONS  /ACRE 


38 

DAT  A 

7?  7 

30 

DATA 

/■  c -i  r 

40 

DAT  A 

cr  h -~i  cr 
? *T  !<  O . 

50 

DAT  A 

4 . 5 ? y . 5 

60 

DAT  A 

7 - 7 !•  * 7 ? cl 

70 

DATA 

•3!  .1  cl  u 7 ? d. 

80 

DAT  A 

3 . 5 ? £ ; 1 

90 

DATA 

y j 1.5’  l 

READY 

RUN 


4 ? 4 j y « 5 ? y 

y ? y . y ? y . 5 ? y 

OJC.-JJC.OJCJ  i ■ •-» 

u?  y?  1 . y?  i . 5 ? 1 
i 5 . y ? . 7 j .5 
. y j . y ? . 4 ? . y 
.5) . 4 ? . y ? . y 


H0M0G1  i£:4£  RTL  THE  06.-' 1 S--"74 


7 SAMPLES 

WITH  7 OBSERUAT I ON 

S EACH. 

SAMPLE 

SUM  OF  SQUARES 

MEAN  SQUARE 

LOG  OF  MEAN 

1 

nr-  »Z|OOcr“7 

. yyyu95 

-b. sol IsE-y 

C. 

6 . W 77 1 4- 

1.01 386 

5. 5483E-3 

S' 

6 . 077 1 4 

1.01 £86 

5. 5483E-3 

*4 

*4  a 2>4o5  r 

. 75S095 

. 1 303  ft. 

cr 

. 1 

7 . 00 

. 98 

— y . 77  y 9 y E. —y 

t« 

4. 34357 

. 708095 

- . 1 49903 

r 

c! > 7 If  14 

. 469534 

a 7 c!  *i'  7 4 cl 

MEAN  OF  ST 

£ - .347075  LOG 

DF  MEAN  St£  = 

~ i a ciy  f 7 ciEl~c. 

DIFFERENCE  = 9.68581E-3 
UNCORRECTED  CHI  SQUARE  = 1.33815 
DEGREES  OF  FREEDOM  = 6 

CORRECTED  CHI  SQUARE  = 1 . 35838 


RUNNING  TIME:  1.1  SECS  L-0  TIME  : 

READY 
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SQUARE 


EXHIBIT  6 


y = production  rate 
x i = slope 
x2  = fuel  loading 
x 3 = (slope) 

Test  H0:  B2=  0 at  95%  CL 


2 

X4  = (fuel  loading) 
x 5 = (slope)  (fuel  loading) 


also  test  H0:  Bi=  B2=  B3=  B4=  B5=  0 


ANQVA 

Source  SS  df 

Regression  due  to  Xi,  x2,  x3,  x4 , X5  857.64370  5 

Regression  due  to  x3 , x3,  x4,  X5  855.57357  4 

Gain  due  to  x2  after  x3,  x3,  x4,  X5  2.07013  1 

Residual  278.5663  141 

Total  1136.21  146 


MS  F 

171.52874  86.822 

2.07013  1.04783 

1 .97563 


F . 0 5 , 1/141  = 3.91  - Accept  H0:  B2  = 0,  F.05,  5/141  = 2.28, 
Reject  H0:  B 1 = B2  = B3  = B4  = B5  = 0 


Test  H0:  B3  = 
Also  test  H0: 

0 at  95%  CL 

Bi  = B3  = B4  = B5  = 

after  eliminating  x2 
0 

AN0VA 

Source 

SS  df 

MS 

F 

Regression  due 
Regression  due 
Gain  due  to  x3 
Residual 
Total 

to  Xi,  x 3 , x4 , X 5 
to  Xi , X4  , X 5 

after  x3 , x4 , x5 

855.57357  4 

847.6201  3 

7.95347  1 

280.63643  142 

1136.21  146 

213.89339 

7.95347 

1.97630 

108.229 

4.02442 

F.o5,  1/142  = 

3.91  - Reject  H0 : 

b3 

= 0 

F.o5,  4/142  = 

2.43  - Reject  H0 : 

Bi 

b3  = b4  = b5 

0 

Test  H0:  B4  = 0 at  95%  CL 

after  eliminating  x2 
AN0VA 

Source 

SS 

df 

MS  F 

Regression  due  to  x1?  x3,  x4,  x5 

855.573 

4 

Regression  due  to  xx,  x3,  x5 

650.293 

3 

Gain  due  to  xu  after  x, , x3,  xs 

205.280 

1 

205.28  103.871 

Residual 

280.63643 

142 

1 .97630 

Total 

1136.21 

146 

F.05,  1/142  = 3.91  Reject  H0:  B4  = 0 
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EXHIBIT  6 (CONTINUED) 


Test  H0 : B x = 0 at  95%  CL  after  eliminating  x2 


ANOVA 


Source 

SS 

df 

MS  F 

Regression  due  to  x19  x3,  x4,  x5 

355.57357 

4 

Regression  due  to  x3,  x4,  x5 

743.712 

1 

Gain  due  to  x,  after  x3,  x4,  xs 
Residual 

111.862 

1 

111.862  56.602 

280.63643 

142 

1 .97630 

Total 

F.05,  1/142  = 3.91  Reject 

1136.21 
H0:  Bj  = 0 

146 

Test  H0:  B5  = 

0 at  95%  CL 

after  eliminating  x2 
ANOVA 

Source 

SS 

df 

MS  F 

Regression  due 

to  X j j X 3 5 x 4 3 x 3 

855.57357 

4 

Regression  due 

to  Xj,  X 3 , X4 

826.672 

3 

Gain  due  to  x5 

after  x1 , x3 , x4 

28.902 

1 

28.902  14.624 

Residual 

280.63643 

142 

1.97630 

Total 

1136.21 

146 

F.o5,  1/142  = 3.91  Reject  H0:  B5  = 0 
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EXHIBIT  7 


READY 

LORD  whwSTEPWI 
READY 

MERGE  smsssTEPW  I ? DATA 
READY 

i 083  D I N Z(10j  1 0 3 j A(i  9 3 ,i  L ( 1 0 3 <■  E ( 1 0 3 ? p t 1 0 3 ? v ( 1 30 , 1 0 3 

i i 40  I F M > 1 30  T HEN  ■ 3 700 

RUN 


SffiSSTLh'WI 


09 ; 39 


0001S074 


TYPE  THE  NUMBER  OF  THE  COLUMN  CORRESPONDING  TO  THE 
DEPENDENT  UAR I ABLE ( Y 3 . ? 1 

DO  YOU  WISH  TO  OMIT  A UAR I ABLE C WK 3 FROM  THE  ANALYSIS?© 
STEP  1 

UAR I ABLE  SELECTED  IS  ...  M 6^ 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP 
F'RuPuRI  lull  UF  UAR  I ABLE  UF  Y REdUL-EB 

PARTIAL  F (D.F.  = 1?  145  3 

CUMULATIUE  SUM  OF  SQUARES  REDUCED.  . 


CUMULATIUE  PROPORTION  REDUCED. .... 
MULTIPLE  CORELATION  COEFFICIENT. . . 

F FuF:  ANALYSIS  UF  UAR.  L D.F.  — 1 ? 1 
STANDARD  ERROR  OF  ESI  I MATE. .... 

UAR I ABLE  REG.  COEFF.  STD. ERR-COEFF 

6 —4 . U4S6 1 E~ 3 3 . 3 4 1 U7E~ 4 

INTERCEPT (A 3 6.66894 


...  .518335 
. . . . 1 56 . 039 
...  5SS.937 

...  .518335  (OF  1136.21 
...  .71 9955 
45  3 156.039 

m m a i ■ !•'  C..  I* 

COMPUTED  T 
- 1 2 . 49 1 6 


*3 


£ ^ V 


. 53.4031 
■ *“r . r 6^0 1 c E c 
..  15.5711 
\ THE  REGRE S S I U N ' 1 


. 565 
.751 


j>S7  (OF  1136.21  3 
989 


1-4  3 93 . 6454 
. 1.85193 


UAR  I ABLE  SELECTED  IS  ...  K 2 
SUM  OF  SQUARES  REDUCED  IN  THIS  STEP. 

PROPORTION  OF  UAR I ABLE  OF  Y REDUCED. 

PART  I FiL  F 1.  D.  F.  — 1 ? 144  J .........  . 

DO  YOU  WISH  TO  ENTER  THIS  UAR I ABLE  I 
CUMULATIUE  SUM  OF  SQUARES  REDUCED. . . 

CUMULATIUE  PROPORTION  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT 

F FOR  ANALYSIS  OF  UAR.  (D.F.  = 2 1 

STANDARD  ERROR  OF  ESTIMATE 

UAR I ABLE  REG.  COEFF.  STB. ERR-COEFF.  COMPUTED  T 

6 •’  — 8.9S7U3E— S - 4.1 98 IE.— 4 — 6..  8 7887 

8 — 7.34U38E— 8 1.S6U18L— 8 —8 . S46U8 

INTERCEPT (A3  7.43566 

STEP  3 & 

UAR I ABLE  SELECTED  IS  ...  X 5 ^ > 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP....  205.28 
PROPORTION  OF  UAR I ABLE  OF  Y REDUCED....  .130671 

F'HF:T  IHL  F (D.F.  — 1 ’ i43  J lUl . 719 

DO  YOU  WISH  TO  ENTER  THIS  UAR I ABLE  IN  THE  REGRESSION?! 

CUMULATIUE  SUM  OF  SQUARES  REDUCED 847.62 

CUMULATIUE  PROPORTION  REDUCED 746008  (OF  1136.21  3 

MULT  I PLE  CORELAT I ON  COEFF  1 0 1 ENT . 8637 1 7 

f for  analysis  uf  uar.  i. d.f.  — s ? 143  j 140.003 
STANDARD  ERROR  OF  ESTIMATE 1.4206 


UAR I ABLE 

6 


INTERCEPT (A3 


REG.  COEFF.  STD. ERR-COEFF. 

2 . 3 1 329E-3  6 . 1 1 283E-4 
~. 83061s  8. 1 1 S25E-8 
— 8 . 8 4 8 1 6E— 3 8 . 8 8 90 8 E— 4 
1 O . 0585 


COMPUTED  T 
3 . 78429 
-10.91 27 
- 1 0 . 0856 
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' 


:■  a 

THIS  STEP. 
REDUCED. 


STEP  4 , .,^C\ 

UAR I flSLE  SELECTED  IS  ... 

SUM  OP  SQUARES  REDUCED  IM 
PROPORTION  OF  URRIABLE  OF  V 
PARTIAL  F i. D.  F . = 1 ? 1 4c!  J ..........  . 

DO  YOU  WISH  TO  ENTER  THIS  OAF:  I ABLE  IN 
CUMULATIUE  SUM  OF  SQUARES  REDUCED......  855.573 

CUMULATIVE  PROPORTION  REDUCED .753808  (OF 

MULTIPLE  CORELATION  COEFFICIENT........  .86776 

f fur  Final Ys I s of  oar.  i. d.f.  — 4 ? 142  j 1 us.  els 


. 7.95347 
. 7' . 0808 IE-3 
. . 4 . 02448 
THE  REGRESSION?! 


1136.E1  J 


STANDARD  ERROR 
VARIABLE 

6 


OF  ESTIMATE i. 

REG . COEFF . STD . ERR-CGEFF . 

£ . 3 i 329E-3  6 . 84924E-4 

— . £98 1 4 1 3 . 96£S7E— £ 

-•£ . £4S  1 6E~3  £ . £U5S7E.— 4 

1 . 9U7U5E— 3 3= 38685E— 4 

18.4141 


1858 1 

COMPUTED  T 
3. 8£4l 


-10. 191 
£. 0861 


{wjA  * 


INTERCEPT (A) 

STEP  5 

OAR I ABLE  SELECTED  IS 

SUM  OF  SQUARES  REDUCED  IM  THIS  STEP... 
PROPORTION  OF  OAF:  I ABLE  OF  Y REDUCED... 
PARTIAL  F (i>.  F . = 1 ? 141  J ...........  . 

DO  YOU  WISH  TO  ENTER  THIS  OAR I ABLE  IN 
CUMULATIUE  SUM  OF  SQUARES  REDUCED..... 

CUMULAT I UE  PROPORT  I ON  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT 

F FOR  ANALYSIS  OF  OAR.  (D.F.  - 5 ? 141 

STANDARD  ERROR  OF  ESTIMATE 1. 

OAR I ABLE  REG.  COEFF.  STD. ERR-COEFF. 


078 1 3 
82197E-3 
. . 1.04  i'33 
THE  REGRESSION?! 
. 857.643 
. .75483  (OF 
. . 868809 
J S6.S££ 


1136.21  j 


4 Hr 


COMPUTED 


;i 1 


304£0S 
i . o60 7 1 L~3 
. . 90705E-3 
04A4Fi;3 


INTERCEPT ( A J 


18. 


593c.  r 
70 1 1 7 
08643 
82364 


'HE 


b . -jbcf'  oc.~4 

6. l48lL— 4 — 

9 . 3 0 46 5 E — 4 
53  3. 95189E-2 

9424 

DO  YOU  WISH  TO  PRINT  THE  TABLE  OF  RESIDUALS?© 

DO  YOU  WISH  TO  COMF'U  i E MORE  REGRESSION?! 

TYPE  THE  NUMBER  OF  THE  COLUMN  CORRESPONDING  TO 
DEPENDENT  OAR I ABLE (Y) . ?1 

Du  YuU  WISH  IU  UMiT  A OAR  I ABLE  (XkJ  FROM  THE  ANALYSIS’:'! 
WHEN  REQUESTED?  ENTER  THE  COLUMN  CORRESPONDING  TO  THE 
OAR I ABLE  TO  BE  OMITTED.  TYPE  O TO  TERMINATE  THE  REQUES 
OAR  I ABLE?  PLEASE’!® 

OAR I ABLE?  PLEASE?© 

STEP  1 

OAR I ABLE  SELECTED  IS  ...  X 6 
SUM  OF  SQUARES  REDUCED  IN 
PROPORTION  OF  OAR I ABLE  OF 
PARTIAL  F (D.F.  = 1?  145  j 

CUMULATIUE  SUM  OF  SQUARES  REDUCED 588.931 

CUMULATIUE  PROPORTION  REDUCED 51333! 

MULT I RLE  CORELAT I ON  COEFF I C I ENT .71 9955 

F F UR  ANALYSIS  OF  UAR.  (D.F.  = 1 ? 145  J 156.839 

STANDARD  ERROR  OF  ESTIMATE. ...... 1 . 94275 

UAR I ABLE  REG.  COEFF.  STD. ERR-COEFF.  COMPUTED  T 

6 —4 „ 0486 1 E~3  3 . £4 1 U7E-4  — 1 £ . 49 1 6 

INTERCEPT (A)  6.66894 

TYPE  THE  NUMBER  OF  THE  COLUMN  CORRESPONDING  TO  THE 
DEPENDENT  UAR I ABLE (YJ . ?1 

DO  YOU  WISH  TO  OMIT  A UARIABLElXKJ  FRUM  I HE  ANFlL  1 Si o : i 
WHEN "REQUESTED?  ENTER  THE  COLUMN  CORRESPONDING  TO  THE 
UAR I ABLE  TO  BE  OMITTED.  TYPE  8 TO  TERMINATE  THE  REQUES 
UAR  I FiBLE  ? PLEASEY’3 
UAR I ABLE  ? PLEASE?© 


TS. 


THIS  STEP. 
Y REDUCED. 


«r.  nr 

•J  1 OOOJ 

1 56 - 039 


(UF  1136.21  J 


TS. 
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STEP  £ 

UFlR I ABLE  SELECTED  IS  ...  X £ 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP....  53.4031 

PROPORTION  OF  UAR I ABLE  OF  V REDUCED 4.7O01£E-£ 

PARTIAL  F (D.F.  = 1?  144  j ......  15.5711 

DO  YOU  WISH  TO  ENTER  THIS  OAR I ABLE  IN  THE  REGRESSION?! 
CUMULAT I UE  SUM  OF  SQUARES  REDUCED......  64£.34 

CUMULATIUE  PROPORTION  REDUCED.. .565337  (OF  H36.£l  ) 

MULTIPLE  CORELATION  COEFFICIENT........  .751889 

F FOR  ANALYSIS  UF  UAR.  i. D.F.  = £ ? 144  j 93.6454 
STANDARD  ERROR  OF  ESTIMATE. ............  1 . 85193 

UAR I ABLE  REG.  CQEFF.  STD. ERR-COEFF . COMPUTED  T 

6 ~ £ . 9 £ 7 U S E ~ 3 4 . 1 9 S 1 E ~ 4 —6. 9 7 £ £ 7 

I NTERCEPT C A3  7 . 43566 
STEP  3 

UAR I ABLE  SELECTED  IS  ...  X 5 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP £85. £8 

PROPORTION  OF  UAR I ABLE  OF  Y REDUCED . 18G671 


PARTIAL  F (D.F. 


1 > 1 43  3 


101, 


19 


DO  YOU  WISH  TO  ENTER  THIS  UAR I ABLE  IN  THE  REGRESSION?! 
CUMULATIUE  SUM  OF  SQUARES  REDUCED......  847. 6£ 

CUMULATIUE  PROPORTION  REDUCED. .........  . 746608  (OF  1136 

MULTIPLE  CORELATION  COEFFICIENT. 863717 

F FOR  ANALYSIS  OF  UAR.  (D.F.  = 3 ? 143  3 140.003 


:i  3 


STANDARD  ERROR  OF  ESTIMATE. 
UAR I ABLE  REG.  COEFF. 

6 £.3l3£9L— 3 

£ -.£30613 

5 — £ . £ 4 S 1 6 E — 3 

10.O585 


STD. ERR-COEFF 


6.  11 
£.  1 1 


388E-4 
3£5E-£ 
-imsE — 4 


4£06 

COMPUTED  T 

-10. 91 £7 
- 1 0 . 0856 


INTERCEPT (A3 
STEP  4 

UAR I ABLE  SELECTED  IS  ...  X 4 
SUM  OF  SQUARES  REDUCED  IN  THIS  STEP. 
PROPORTION  OF  UAR I ABLE  OF  Y REDUCED. 
PARTIAL  F (D.F.  = 1?  14£  3. 


DO  YOU  WISH  TO  ENTER  THIS  UAR I ABLE 
CUMULATIUE  SUM  OF  SQUARES  REDUCED. 

CUMULATIUE  PROPORTION  REDUCED 

MULTIPLE  COREL AT  I ON  COEFFICIENT. . . 
F FOR  ANALYSIS  OF  UAR.  (D.F.  - 4 ? 

STANDARD  ERROR  OF  ESTIMATE........ 

UAR I ABLE  REG.  COEFF. 

6 „ £.313£9E-3 

£ -.£98141 

5 — £ . £ 4 816 E — 3 


T m 

1 I I 


. 7.95347 
. 7.G0001E-3 
. . 4 . UE44E 
THE  REGRESSION?! 


(OF  1 136. £1 


1 4£ 


4 


1 . 90705E-3 


STD. ERR-COEFF 
6. U49£4E-4 
8 . 9886  E £ 
£ . £0587E-4 
9 . 5U8E5E-4 


. i'  530tiy 
. 8677 6 
3 108. ££9 
1 . 4058 1 

COMPUTED  T 
3 . 8E4 1 


i.  1917 
006 1 


INTERCEPT (A3  10.4141 

DO  YOU  WISH  TO  PRINT  THE  TABLE  OF  RESIDUALS?© 
Du  YuU  W 1 8 FI  To  OUMF'UTE  MURE  REGRESSION’:' 1 
STEP  1 

UAR I ABLE  SELECTED  IS  ...  X 6 
SUM  OF  SQUARES  REDUCED  IN  THIS  STEP 
PROPORTION  OF  UAR I ABLE  OF  Y REDUCED 
PARTIAL  F (D. F.  = 1 ? 145  3 ........  . 

CUMULATIUE  SUM  OF  SQUARES  REDUCED.. 

CUMULAT I UE  PROPORT  I ON  REDUCED 

MULT  I RLE  CQRELAT I ON  COEFF I C I ENT . 

F FOR  ANALYSIS  OF  UAR.  (D.F.  = 1 
STANDARD  ERROR  OF  ESTIMATE. 

UAR I ABLE  REG.  COEFF. 

6 -4 . 0486 1 E-3 

I NTERCEPT ( A 3 6 . 66894 


1 


STD. ERR-COEFF 
3 . £4 1 07E-4 


800. 

. 5 1 8335 
. 1 36 . 039 

588.937 
. 5 1 8335  ( OF 
. 7 1 9955 
3 156.039 

i ‘~T » 4- 1~‘  cj 

COMPUTED  T 
- 1 £ . 49 1 6 


1 1 36 . £ 1 3 


-41- 


mm 


1. 


mam 


STEP  £ 

UAR I ABLE  SELECTED  IS  ...  X £ 

SUM  OF  SQUARES  REDUCED  III  THIS  STEP....  53.4631 
0 

NOW  FIT  3416 


ready 

GO  TO  1 

TYPE  THE  NUMBER  OP  THE  COLUMN  CORRESPONDING  TO  THE 
DEPENDENT  UflR  I ABLE  f Y ) . ? 1 

DO  YOU  WISH  ro  OMIT  Fi  UARI ABLE  (XKJ  PF'OM  THE  ANALYSIS?! 

WHEN  REQUESTED?  ENTER  THE  COLUMN  CORRESPONDING  TO  THE 

UARI ABLE  TO  BE  OMITTED.  TYPE  6 TO  TERMINATE  THE  REQUESTS. 

UARI  ABLE?  PLEASE  YlN 

UARI  ABLE?  PLEASE/\5y 

UARI ABLE?  PLEASE ?0 

STEP  1 

UARI ABLE  SELECTED  IS  ...  X 6 
SUM  OF  SQUARES  REDUCED  IN  THIS  STEP. 

PROPORTION  OF  UARI ABLE  OP  Y REDUCED 

PARTIAL  F (D.F.  = 1?  145  3 .... 

CUMULAT I UE  SUM  OF  SQUARES  REDUCED. . 

CUMULAT I UE  PROPORT I ON  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT. . . . 

F FuF!  ANALYSIS*  OF  UflR.  UJ.F.  -■  i ? 

STANDARD  ERROR  OF  ESTIMATE... 

UARI ABLE  REG.  COEFF.  STD. ERR-COEFF. 

6 —4. U4861L— 3 . 3.  £41U7E~4 

I NTERCEPT ( A 3 6 . 66894 

S I E.P  £ 

UARI ABLE  SELECTED  IS  ...  X £ 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP 53.4031 


5 i 0335  top  1136, 
7 1 9955 
145  3 156.039 
. . . . 1 . 94£75 

COMPUTED  T 
~1£. 4916 


REDUCED. 


PROPORTION  OF  UARI ABLE  OF  V 
PARTIAL  F i. D.P.  — 1?  144  3 . 

DO  YOU  WISH  TO  ENTER  THIS  UARI ABLE  If 
CUMULAT I UE  SUM  OF  SQUARES  REDUCED. 

CUMULAT I UE  PROPORTION  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT... 
f f uF:  final Ys  i s of  uar  . i D . f . — £ ? 

STANDARD  ERROR  OF  ESTIMATE .... 

UAR I ABLE  REG.  COEFF.  STD. ERR-COEFF 

S — £ . S1  £7  U 3 E — 3 4 . 1 9 £ 1 E — 4 

£ —? . 34U3£E~£  1 . £.60 1 £E~ £ 


. 4 . 700 1 £E-£ 

..  15.5711 
THE  REGRESSION?! 
. 64 £. 34 


5 £53 3 f 
751889 


l i 


136. £1  3 


4 3 


1.85193 

COMPUTED  T 

-3. 9460£ 


I NT  ERuEF'T  l A J 4 7 . 43566 
STEF'  3 

UAR I ABLE  SELECTED  IS  ...  X 4 
SUM  OF  SQUARES  REDUCED  IN  THIS  STEP. 
PROPORTION  OP  UAR I ABLE  OF  Y REDUCED. 

PART  I FlL  F 1.  D.  F.  = 1 ? 143  .1  .............  . £.  3* 

DO  YOU  WISH  TO  ENTER  THIS . UARI ABLE  IN  THE  REC 

CUMULAT I UE  SUM  OP  SQUARES  REDUCED 650.E 

CUMULAT I UE  PROPORT I ON  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT.. 

F FOR  ANALYSIS  OF  UAR.  ( D.F.  = 3 
STANDARD  ERROR  OF  ESTIMATE....... 

UAF:  I ABLE.  REG . CoEF  F . 

6 ~£»  9£703E-3 


143 


. fjK 

3 63.  1 
1 . 84S 
STD. ERR-COEFF.  ( 
4.1  f'£6yE-4 


95347 
0000 1 E-3 
■ '063 

RESS I ON? 1 
93 

37  UJF  11 36. 
£9 
917 


:1  3 


4 

INTERCEPT  I A 3 


- . 1 4093 1 
1 . 96765E-3 
r . 79 1 £8 


. 04  1"  fclfciD 

1. £4651 E-3 


1 

OMPUTED  ' 
-7. 00465 
-£. 94435 
1 . 5£99 1 
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' 


K 


THIS  STEP 
Y REDUCED 


■-T- 


39 


145  ) 


O i boo 

1 56 . 8 

•So  c So? 

518335 
719955 
1 56 . 039 
» 94475 


(OF  1136.41  J 


no  YOU  WISH  TO  PRINT  THE  TfiELE  OF  RESIDUALS?© 

DO  YOU  WISH  TO  COMPUTE  MORE  REGRESSION?! 

TYRE  THE  NUMBER  OF  THE  COLUMN  CORRESPONDING  TO  THE 
DEPENDENT  OAR I ABLE ( Y ) . ? 1 

DO  YOU  WISH  TO  OMIT  A URR I ABLE ( ?K ) FROM  THE  ANALYSIS?! 
WHEN  REQUESTED.-  ENTER  THE  COLUMN  CORRESPONDING  TO  THE 
WAR I ABLE  TO  BE  OMITTED.  TYPE  0 TO  TERMINATE  THE  REQUESTS 
WAR I ABLE-  PLEAS 
WAR I ABLE-  PLEAS 
OAR I ABLE-  PLEAS 
STEF'  1 

UARIABLE  SELECTED  IS  ...  > 

SUM  OF  SQUARES  REDUCED  IN 
F'PUF'URT  I UN  UF  UARIABLE  UF 
F'ARTIAL  F i. D.  F . — 1 - 145  J ......  . 

CUMULAT I UE  SUM  OF  SQUARES  REDUCED 

CUMULAT I UE  PROPORT I ON  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT.. 

F FOR  ANALYSIS  OF  UAR . (D.F.  = 1 

STANDARD  ERROR  OF  ESTIMATE... 

UAR I ABLE  REG . CUEFF 

6 -4. 0486 1E- 

INTeRcEPT  i.  A J 6.  66 6 94 
' STEP  4 

UARIABLE  SELECTED  IS 
SUM  UF  SQUARES  REDUC 
PROPORTION  OF  UARIABLE  OF 
PART  IAL  F i.  D.  F . = 1 - 144  1.  ....... 

DO  YOU  WISH  TO  ENTER  THIS  UARIABLE 
CUMULAT I UE  SUM  OF  SQUARES  REDUCED. 

CUMULAT I UE  PROPORTION  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT... 

F FOR  ANALYSIS  OF  UAR.  (D.F.  = 4 - 
ERROR  OF 
RE( 


STD. ERR-COEFF. 
3. E4107E-4 


COMPUTED  "I 
-14. 4916 


ED  IN 


THIS  STEP. 
r‘  REDUCED. 


IN 


•1 


S rAMDARD 
UARIABLE 


. 31.8985 
. 4 . S07 46E-4 
..  8.91475 
THE  REGRESSION" 

....  bet! . bb-J 

....  . 5464 1 i.  UF  1 1 66 
....  . 7 69 1 95 
144  j 86.7337 

ESTIMATE 1.89184 

. COEFF.  STD. ERR-COEFF.  COMPUTED  T 


£1  J 


INTERCEPT (A) 

STEP  3 
UARIABLE  SELECTED  IS 
SUM  OF  SQUARES  REDUC 


-6.  64966L— 6 


6 . 93968 


1 3883E-4 
i 94867E-4 


- r . 85069 
-6.98544 


ED 


IN  THIS  STEF'.. 
PROPORTION  OF  UARIABLE  OF  Y REDUCED.. 

PARTIAL  F (D.F.  = 1-  143  j 

DO  YOU  WISPUTO  ENTER  THIS  UARIABLE  IN 
CUMULAT I UE  SUM  OF  SQUARES  REDUCED 

CUMULAT I UE  PROPORTION  REDUCED 

MULTIPLE  CORELATION  COEFFICIENT.. 

F FOR  ANALYSIS  OF  UAR.  (D.F.  = 3 

STANDARD  ERROR  OF  ESTIMATE 

UARIABLE 


143 


REG.  COEFF.  STD. ERR-COEFF. 


■ 166.6  !*  r" 

. .108147 
..  44.7685 
THE  REGRESSION?! 

. f'  4 6 • r 1 c’ 

. .654557  (OF  113 
. . 8O9047 
90 . 3£03 
1 . 65674 

COMPUTED  T 


6.41  J 


INTERCEPT (A) 


6.69155E— . 
—4. 16795E— 
— 1 . 606 1 9E— 4 
6.  1 666-5 


6. 66747E-4 
5. 91 193E-4 
6 . S9456E— 4 


.361587 
-7 . 05O05 
-6. 69093 
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DO  YOU  WISH  TO  PRINT  THE  TflBLE  OF  RES  I BURLS'' O 
DO  YOU  WISH  TO  COMPUTE  MORE  REGRESSION' 1 
TYPE  THE  NUMBER  OF  THE  COLUMN  CORRESPONDING  TO  THE 
DEPENDENT  UfiR I RELE ( Y ) . ? 1 

DU  YuU  WISH  TU  UNIT  H UFlRIRELE  i.WRj  FROM  THE.  ANALYSIS? 1 
WHEN  REQUESTED ? ENTER  THE  COLUMN  CORRESPONDING  TO  THE 
UFlRIRELE  TO  EE  OMITTED.  TYPE  O TO  TERMINRTE  THE  REQUESTS. 


UfiR  I FlELE  ? PLEfiSE93\ 

UfiR  I RBLE  ? PLEASE  ^R/ 

UFlRIRELE?  PLEASE' O 
s i EP  i 

UFlRIRELE  SELECTED  IS  ...  X,  £ 

SUM  OF  SQURRES  REDUCED  IN  THIS  STEP....  475. Sib 

PROPORTION  OF  UfiR I RELE  OF  Y REDUCED .4186 

F'RFTIRL  F l D.  F ■ = 1 ? i 40  J .............  . 1W4.  398 

CUMULRT I UE  SUM  OF  SQUARES  REDUCED......  475. Si 6 

CUMULRTIUE  PROPORTION  REDUCED 4186  (OF  1136.81  ) 

MULTIPLE  CORELRT I ON  COEFFICIENT .646993 


F FOR  ANALYSIS  OF  URR.  (D.F.  = 1 ? 145  j 104.398 

s i rndrrd  errur  uf  estimate. ............  e. is44s 

URR I ABLE  PEG.  COEFF.  STD. ERR-COEFF.  COMPUTED  T 

£ - . 1 6 1 £ 1 4 1 . 5778£E-£  -18. £175 

INTERCEPT (R)  7.43566 

STEP  £ 

URR I RELE  SELECTED  IS  ...  X 5 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP 343. 10£ 


PROPORTION  OF  UFlRIRELE  OF  Y REDUCED. 

PARTIAL  F i.  D.F.  = 1?  144  .i  . . 

DO  YOU  WISH  TO  ENTER  THIS  UfiR I RELE 
CUMULRTIUE  SUM  OF  SQURRES  REDUCED 
CUMULRTIUE  PROPORTION  REDUCED.  . . . 
MULTIPLE  CORELRT I ON  COEFFICIENT.. 

F FOR  ANALYSIS  OF  URR.  (D.F.  = £ 

s i hndrrd  errofi  uf  estimate.  ...... 

UFlRIRELE  REG.  COEFF. 

£ - . 1 6 1 £ 1 4 

5 -1.53116E-3 

INTERCEPT (R)  9.£££01 

STEP  3 

URR  I FlELE  SELECTED  IS  ...  X 4 
SUM  OF  SQURRES  REDUC 


IN 


144 


STD. ERR-COEFF 
1 . 89763E-E 
l.££74£E-4 


. SO 1 S7£ 

. 155.617 

HE  REGRESSION?! 

818. 719 

. 78857 1 i.  UF  1 1 S6 . 8 1 
. 848865 
j 185.669 
1 . 48485 
. COMPUTED  T 
— 1 4. oS7 4 
id.  a 4 !*'  4 f' 


ED 


THIS  STEP. 
Y REDUCED. 


IN 

PR  OP  OF:  I ION  OF  URR  I ABLE  OF 

F'RFlTIRL  F i.  D.  FK  = 1 ? 1 4S  .i  .....  

DO  YOU  WISH  TO  ENTER  THIS  URR I RBLE  IN 
CUMULRTIUE  SUM  OF  SQURRES  REDUCED.  . . . 

CUMULRT I UE  PROPORT I ON  REDUCED 

MULTIPLE  CORELRT  I ON  COEFFICIENT 

F FOR  ANALYSIS  OF  URR.  (D.F.  = 3 ? 143  ) 


. 7.95347 
. 7. 8800 IE -3 
. . 3 . 67 436 
THE  REGRESSION'"  1 
. 886 . 678 

. . 7£757 1 ( OF  1 1 36 , 


ii  J 


8 -i  89  r'  8 

1 £7 . 383 


STANDARD 
URR I RELE 


4 


INTERCEPT (Rj 


ERROR  OF 

REG.  COEFF. 
— . £887 48 
-1.53116E-3 
1 . 9 87 USE -3 
=■ . 57763 


ESTIMATE 1.471  £5 

STD. ERR-COEFF.  COMPUTED  T 
3 . 6869 1 E— £ —6 . £04 i 7 

1 . E1617E-4  -18.59 

9 . 94879E— 4 1 . 9 1 686 
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DO  YOU  WISH  TO  PRINT  THE  TABLE  OF  RESIDUALS?© 

DO  YOU  WISH  TO  COMPUTE  MORE  REGRESSION?! 

TYPE  THE  NUMBER  OF  THE  COLUMN  CORRESPONDING  TO  THE 
DEPENDENT  UAR I ABLE ( Y) . ? 1 

DO  YOU  WISH  TO  OMIT  A UAR I ABLE (XK)  FROM  THE  ANALYSIS?! 

WHEN  REQUESTED ? ENTER  THE  COLUMN  CORRESPONDING  TO  THE: 

UAR I ABLE  TO  BE  OMITTED.  TYPE  8 TO  TERMINATE  THE  REQUESTS. 

UAR  I ABLE?  PLEASE'TSj 
UAR  I ABLE?  PLEASEW 
UAR I ABLE?  PLEASE78 
STEP  i 

UAR I ABLE  SELECTED  IS  ...  X 6 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP 588.937 

PROPORTION  OF  UAR I ABLE  OF  Y REDUCED 518335 

PARTIAL  F (D.F.  = 1?  145  ) . . 156.839 

CUMULATIUE  SUM  OF  SQUARES  REDUCED......  588.93? 

CUMULATIUE  PROPORTION  REDUCED 518335  (OF  1136.21  ) 

MULT  I RLE  CORELAT I ON  COEFF I C I ENT 71 9955 

F F UR  ANALYSIS  UF  UAR.  i. D.F.  — 1 ? 145  i 156.839 

STANDARD  ERROR  OF  ESTIMATE 1.94275 

UAR I ABLE  REG.  COEFF.  STD. ERR-COEFF.  COMPUTED  T 

6 —4 . 8486 1 E— 3 3 .241 87E—4  — 1 c . 49 1 6 

I NT  ERCEF'T  •.  A J 6 . 66894 
STEP  2 

UAR I ABLE  SELECTED  IS  ...  X 2 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEP 53.4831 

PROPORTION  OF  UAR I ABLE  OF  Y REDUCED 4.78812E-2 

PARTIAL  F (D.F.  = 1?  144  ) . 15.5711 

DO  YOU  WISH  TO  ENTER  THIS  UAR I ABLE  IN  THE  REGRESSION?! 

CUMULATIUE  SUM  OF  SQUARES  REDUCED 642.34 

CUMULATIUE  PROPORTION  REDUCED .565337  (OF  1136.21  j 

MULT  I RLE  CORELAT  I ON  COEFF  I C I ENT 751 889 

F FOR  ANALYSIS  OF  UAR.  (D.F.  = 2 ? 144  j 93.6454 

STANDARD  ERROR  OF  ESTIMATE. 1.85193 

UAR I ABLE  REG.  COEFF.  STD. ERR-COEFF.  COMPUTED  T 

6 — 8. S2703E— 3 4. 198 IE- 4 — 6.97887 

2 - 7 . 34832E-2  1 . 868 1 8E-2  -3 . 94682 

INTERCEPT (A)  7.43566 

STEP  3 

UAR I ABLE  SELECTED  IS  ...  X 5 

SUM  OF  SQUARES  REDUCED  IN  THIS  STEF'....  285.28 

PROPORTION  OF  ■VARIABLE  OF  Y REDUCED 188671 

PART  I AL  F i.  D.  F . = 1 ? 143  1..............  181.719 

DO  YOU  WISH  T£  ENTER  THIS  UAR I ABLE  IN  THE  REGRESSION?! 

CUMULATIUE  SUM  OF  SQUARES  REDUCED 847.62 

CUMULATIUE  PROPORTION  REDUCED 746888  (OF  1136.21  j 

MULT  I RLE  CORELAT  I ON  COEFF  I C I ENT 8637 1 7 

F FuFl  ANAL Ys 1 8 UF  UAR.  (Jj.F.  — 3 ? 143  J 148.883 

STANDARD  ERROR  OF  ESTIMATE 1.4286 


UAR I ABLE 

REG.  COEFF. 

STD. ERR-COEFF. 

COMPUTED  T 

f, 

C.  a cl  1 -2'  c!  ‘y  El  i' 

6.11 888E- 4 

a »'  d4d  2* 

CL 

- . 2386 1 3 

CL  a i 1 -2'iIL^jE*.  cl 

-18. 9127 

c* 

INTERCEPT (A) 

—8. 848 1 6E— 3 
1 8 . 8585 

8. 88908E— 4 

- 1 8 . 8856 

DO  YOU  WISH 

TO 

PRINT  THE  TABLE  OF  RESIDUALS 

?0 

DO  YOU  WISH 

TO 

COMPUTE  MORE 

REGRESSION?© 

END  OF  STEF'WI 
NOW  AT  END 
READY 


45- 


' 


